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Abstract 

Recent developments in the physics of ultracold gases provide wide possibilities for reduc- 
ing the dimensionality of space for magnetically or optically trapped atoms. The goal of 
these lectures is to show that regimes of quantum degeneracy in two-dimensional (2D) 
and one-dimensional (ID) trapped gases are drastically different from those in three di- 
mensions and to stimulate an interest in low-dimensional systems. Attention is focused on 
the new physics appearing in currently studied low-dimensional trapped gases and related 
to finite-size and finite-temperature effects. 
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INTRODUCTION 



The subject of low-dimensional quantum gases has a long pre-history. The influence of dimen- 
sionality of the system of bosons on the presence and character of Bose-Einstein condensation 
(BEC) and superfluid phase transition has been a subject of extensive studies in the spatially 
homogeneous case. From a general point of view, the absence of a true condensate in 2D 
and ID at finite temperatures follows from the Bogoliubov k~ 2 theorem and originates from 
long- wave fluctuations of the phase (see, e.g., [1]). This has been expounded by Mermin and 
Wagner [2] and by Hohenberg [3], and formed a basis for later investigations. 

The earlier discussion of low-dimensional Bose gases was mostly academic as there was 
no possible realization of such a system. Fast progress in evaporative and optical cooling of 
trapped atoms and the observation of Bosc-Einstcin condensation (BEC) in trapped clouds 
of alkali atoms [4-6] stimulated a search for non-trivial trapping geometries. Present facilities 
allow one to tightly confine the motion of trapped particles in one (two) direction(s) to zero 
point oscillations. Then, kinematically the gas is 2D (ID), and the difference from purely 
2D (ID) gases is only related to the value of the effective interparticle interaction which now 
depends on the tight confinement. 

Recent experiments have already reached 2D and ID regimes for trapped Bose gases and 
studied some of the quantum degenerate states. These studies bring in new physics originating 
from a finite size of the system, spatial inhomogeneity, and finite temperatures. The present 
lectures cover most important issues in the physics of 2D and ID trapped quantum gases: the 
nature of various quantum degenerate states, the role of interaction between particles, and the 
role of finite-temperature and finite-size effects. 



1 LECTURE 1. BEC IN IDEAL 2D AND ID GASES 



We start with describing a cross-over to the BEC regime in ideal 2D and ID Bose gases with a 
finite number of particles. We will consider an equilibrium gas at temperature T in the grand 
canonical ensemble, where the chemical potential /x is fixed and the number of particles N is 
fluctuating. In the thermodynamic limit (N — ► oo) this is equivalent to the description in the 
canonical ensemble (fixed N and fluctuating fi). 

In any dimension and confining potential the gas is characterized by a set of eigenenergies 
of an individual particle, E u , with the index v labeling quantum numbers of the particle 
eigenstates. The (average) total number of particles N is then related to the temperature and 
chemical potential by the equation 

where N u (z) — l/(expz — 1) are the equilibrium occupation numbers of the eigenstates. We 
now demonstrate how Eq. (f 1 . 1 1> allows one to establish the presence or absence of BEC in 2D 
and ID ideal Bose gases. 
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1.1 Uniform ideal gas 

In an infinitely large uniform gas the particle eigenstates are characterized by the momentum 
k and the eigenenergy Ek — h 2 k 2 /2m, where m is the mass of a particle. Then Ea. Hl.l(l takes 
the form 

with d being the dimension of the system, and £1 the d-dimensional volume. 
In the 2D case the integration in Ea. l|1.2(l is straightforward and we obtain 

/ Lt = Tln[l-exp(-n 2 A|)] <0 , (1.3) 

where At — (2irh 2 /mT) 1 / 2 is the thermal de Broglie wavelength, and n 2 is the 2D density. The 
quantity n 2 A|, is called the degeneracy parameter and in 2D it can be written as n 2 A 2 n — T /Td, 
where Td = 2-7r?i 2 n 2 /m is the temperature of quantum degeneracy. In the limit of a classical 
gas, n 2 A|, <C 1, Eq. i|1.3fl gives the well-known result p = Tln(n 2 A^). For a strongly degenerate 
gas, where n 2 A|, 3> 1, we obtain p = — Texp(— n 2 A^). Unlike in the 3D case, the dependence 
p(T) is analytical and it shows a monotonic increase of the chemical potential with decreasing 
temperature up to T — > 0. In the thermodynamic limit the population of the ground state 
(k = 0) remains microscopic. One thus can say that there is no BEC in a finite-temperature 
ideal uniform 2D Bose gas. 

The situation is similar for an infinite uniform ID Bose gas, where the degeneracy parameter 
is niA T — {T/Td) 1 / 2 and the temperature of quantum degeneracy is given by T d = 27rft 2 n 2 /m, 
with n\ being the ID density. In the classical limit (niAy <C 1), and in the limit of a strongly 
degenerate gas (n-i Aj- 3> 1), Eq. i|1.2|) gives 

/i = Tln(niA T ) , mA T < 1 ; (1.4) 
irT 

M = -7 — 7—^ i niA T >l . (1.5) 
(niA T ) 2 

Again, the chemical potential monotonically decreases with temperature and remains negative 
at any T, which indicates the absence of BEC. 

The absence of BEC in infinitely large uniform 2D and ID Bose gases is a striking difference 
from the 3D case. This difference originates from the energy dependence of the density of states. 
The (energy) density of states is p(E) oc £ , ( d / 2 ^ 1 ', where d is the dimension of the system, 
and for the 3D gas it decreases with E. Therefore, at sufficiently low temperatures it becomes 
impossible to thermally occupy the low energy states while maintaining a constant chemical 
potential or density. As a result, a macroscopic number of particles goes to the ground state 
(k = 0), i.e. one has the phenomenon of BEC. In 2D and ID the density of states does not 
decrease with E and this phenomenon is absent. 

1.2 Ideal gas in a harmonic trap 

For 2D and ID Bose gases in a harmonic confining potential the density of states is p(E) oc 
j^(d-i) anc j situation changes. The population of the ground state (E = 0) is 

No = ( l -FF\ i ' ( L6 ) 

exp(-/x/T) - 1 

For a large but finite number of particles in a trap, Nq can become macroscopic (comparable 
with N) at a small but finite negative p. One then speaks of a cross-over to the BEC regime. 
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We first discuss the BEC cross-over for the 2D Bose gas in a symmetric harmonic confining 
potential V(r) = muj 2 (x 2 + y 2 )/2. In this case the particle energy is E v — hu(n x +n y ), with 
quantum numbers n xi n y being non- negative integers. The density of states is then p(E) = 
EJ (fiui) 2 , and the contribution of low-energy excited states to the total number of particles is 
negligible. Therefore, separating out the population of the ground state, one can replace the 
summation in Eq. i|l.l[l by integration: 



N = Ni 



o 



dE p{E)N\ 



E- p 
T 



Assuming a large population of the ground state, from Ea. (|1.6ll we obtain 
and the population of excited trap states proves to be 



-H/T: 



dEp{E)N< 



E 



T 



T 



l + lnA 
A 



This allows us to write Ea. (|1.7() in the form 



A 



= N a - [ — 

nio 



T \ 1 + lnA 



An 



where 




(1.7) 
1/A « 1, 

(1.8) 

(1.9) 
(1.10) 



For a large number of particles, Eq. (|1.9|) indicates the presence of a sharp cross-over to the 
BEC regime at T w T c . Below T c we omit the last term in Ea. l|1.9|l and obtain the occupation 
of the ground state 



An « A 



(1.11) 



(Nmuj 2 /T c ) and the de Broglie wavelength of particles 
yjh 2 jmT c becomes comparable with the mean interparticle separation ~ {Nmuj 2 /T^ 1 / 2 



Note that at T c the particle density is 
A 



The result of Eq. I|l.ll|) is similar to that in the 3D case and it was first obtained by Bagnato 
and Kleppner [7] . 

Above T c one can omit the first term on the rhs of Eq. (|1.9(l . The width of the cross-over 
region, i.e. the temperature interval where both terms on the rhs of Ea. (|1.9fl are equally 
important, is given by 

AT /In A 

The Bose-condensed fraction of particles, An(T)/A, following from Ea. (|1.9|l . is presented in 
Fig- D for various values of A. For a large A the cross-over region is very narrow and one can 
speak of an ordinary BEC transition in an ideal harmonically trapped 2D gas. 

For the ID Bose gas in a harmonic potential V(z) — mLu 2 z 2 /2, the particle energy is 
E = hioj, with j being a non-negative integer, and the density of states is p(E) = 1/fku. Here 
the integral representation of Ea. p.7|l fails as the integral diverges for \i — > 0. Therefore, we 
should correctly take into account the discrete structure of the lowest energy levels. In the 
limit {— /i, Hlu} -C T we rewrite Ea. (|l.l|) in the form 



rp M 



1 



3=1 J 



E 



i 



(1.13) 
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N Q /N 



N=50 
N=200 
N=10 3 
N=10 6 




T/T 



Fig. 1. The ground state population in a 2D trap versus temperature, calculated from Eq. 



where the number M satisfies the inequalities 1 -C M <C T/hu. The first sum is 



M 



l —nr- = i>{M + 1 - n/tvjj) - ip(l - n/hJ) w -ip(l - n/ftJ) + ln(M - n/hu), (1.14) 

3=1 J W 

where ip is the digamma function. The second sum in Ea. i|1.13[) can be transformed to an 
integral 



y, 1 T f°° dx T hw(M — n/hu) 

^ 1 exp(tkjj/T - fi/T) - 1 fko JtujjM/T exp(ir — n/T) — 1 hw T 



j=M+ 



Finally, since in the limit of \n\ «T the chemical potential is related to the population of the 
ground state as — [i w T/Nq, we reduce Ea. (|1.13|l to the form 



N T \ ( T ^ N T (l T 



(1.15) 



As in the 2D case, we have two regimes, with the border between them at a temperature 

N 



In TV 



(1.16) 



For temperatures below T\d, the first term on the rhs of Eq. l|1.15[) greatly exceeds the second 
one and the ground state population behaves as 



(1.17) 



The cross-over region is determined as the temperature interval where both terms are equally 
important: 

1\D m TV 
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Fig. 2. The ground state population in a ID trap versus temperature, found from Eq. Ijl.15^ . 

In contrast to the 3D and 2D cases, the cross-over temperature is much lower than the degen- 
eracy temperature Td ~ Nfku. The described results have been obtained by Ketterle and van 
Druten [8]. In Fig[5]we present the relative occupation of the ground state Nq(T)/N calculated 
from Eq. Ql.15j) . 

The cross-over region in the ID case is much wider than in 2D. This is not surprising as the 
cross-over is present only due to the discrete structure of the trap levels. The quasiclassical 
calculation does not lead to any sharp cross-over [7]. 

Problem: Describe the cross-over to the BEC regime in an ideal 2D gas of N identical 
bosons in a rectangular box with dimensions L Xl L y . Find the cross-over temperature and the 
width of the cross-over region (E.B. Sonin, 1969). 

2 LECTURE 2. INTERACTIONS AND BEC REGIMES IN 2D TRAPPED GASES 

Interaction between particles drastically changes the picture of Bose-Einstein condensation in 
2D and ID gases. In the uniform 2D gas a true condensate can exist only at T — and, as 
mentioned in the Introduction, its absence at finite temperatures originates from long-wave 
fluctuations of the phase [1] . However, as was pointed out by Kane and Kadanoff [9] and then 
proved by Berezinskii [10], there is a superfluid phase transition at sufficiently low T . Kosterlitz 
and Thouless [11] found that this transition is associated with the formation of bound pairs of 
vortices below the critical temperature 

Tkt = ^n s , (2.1) 

where n s is the superfluid density just below Tkt- This temperature is of the order of or 
smaller than Td/4, where Td — 2irh 2 n/m is the temperature of quantum degeneracy introduced 
in Lecture 1, and the notation n is used in this Lecture for the 2D density. Recent Monte Carlo 
calculations [12] established an exact relation between Tkt and Td for the weakly interacting 
2D Bose gas. 

Early theoretical studies of 2D systems have been reviewed by Popov [1] and have led to 
the conclusion that below the Kostcrlitz-Thoulcss transition temperature the Bose liquid (gas) 
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is characterized by the presence of a quasicondensate, that is a condensate with fluctuating 
phase (see, e.g., [13]). In this case the system can be divided into blocks with a characteristic 
size greatly exceeding the healing length but smaller than the phase coherence length. In each 
block one has a true condensate but the phases of different blocks are not correlated with each 
other. 

The Kosterlitz-Thouless transition has been observed in monolayers of liquid helium [14], 
and recently the observation of this transition has been reported for the 2D gas of spin-polarized 
atomic hydrogen on liquid helium surface [15]. The Kosterlitz-Thouless transition is discussed 
in the lectures of Benoit Doucot. 

In this Lecture we focus our attention on the interaction between particles and on the 
correlation properties of 2D Bose-condensed gases. Recently, the 2D regime was realized for 
trapped atomic Bose-Einstein condensates at MIT [16], LENS [17], Innsbruck [18], and JILA 
[19]. This was done by a tight optical confinement of the particle motion in one direction [16,17], 
or by a fast rotation of the cloud leading to an increase in the size of the sample in two 
directions [19]. We will focus special attention on how the nature of the Bose-condensed state 
is influenced by a finite size of the system. 



2.1 Weakly interacting regime 

We will consider weakly interacting gases with a short-range potential of interaction between 
particles. In this case the total interaction energy is equal to the sum of pair interactions and 
can be written as E- lnt = N 2 g/2Q, where g is the coupling constant for the pair interparticle 
interaction, N is the number of particles, and Q the volume. Accordingly, the interaction 
energy per particle is equal to I = ng, with n being the 2D density. 

Let us now discuss the conditions which are required for the 2D gas to be in the weakly 
interacting regime. The commonly used criterion assumes that the mean interparticle separa- 
tion f greatly exceeds the characteristic radius of interaction between particles, R e . In the 2D 
case we have f rs (27m) -1 ' 2 and thus obtain the inequality 

nR 2 e < 1 . (2.2) 

The weakly interacting regime requires that at interparticle distances of the order of f, the 
wavefunction of particles is not influenced by the interaction between them. Relying on this 
requirement, we develop a physical picture which will be used for finding how the criterion of 
the weakly interacting regime depends on the dimensionality of the system. We consider a box 
of size f, which on average contains one particle. In the limit of T — > 0, the particle kinetic 
energy is K « h 2 /mf 2 . The wavefunction of the particle is not influenced by the interparticle 
interaction if K is much larger than the interaction energy per particle, I — ng. The inequality 
K 3> \I\ immediately gives the criterion of the weakly interacting regime in terms of the density 
and coupling constant. In the 2D case this criterion takes the form 

. (,3) 

In the 2D gas at T — ► the coupling constant is (see, e.g. [1]) 

4irh 2 1 



m ln(l/ncZ 2 ) 



(2.4) 



where a length depends on the shape of the interatomic potential and in the absence of 
scattering resonances is of the order of R e . Then, from Eq.J22|) we immediately arrive at the 
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criterion (|2.2|l . However, on approach to a resonance, the value of is quite different from 
R e , and one should return to the criterion l|2.3[l . 

In the dilute limit Ea. (|2.4|) gives g > 0, except for the case where a weakly bound state of 
colliding atoms is present. Then the length d* is extremely large and even at very low densities 
one can have the condition nd 2 3> 1 leading to an attractive mean- field interaction (g < 0). 
However, in this case the 2D Bose gas at T — > is unstable with regard to collapse. 

Note that in the 3D case we have the coupling constant g^jj = ATrh 2 a/m and r ~ n 3 ^ 3 , 
with n 3 D being the 3D density, and a the 3D scattering length. The condition K ^> |/| then 
leads to the well-known criterion n\a\ <C 1, required for the weakly interacting regime in the 
3D gas. In the absence of resonances, the 3D scattering length is of the order of R e and this 
criterion is equivalent to the inequality R e <C f. 

In ongoing experiments, two-dimensional atomic gases are obtained by (tightly) confining 
the motion of particles in one direction to zero point oscillations [16-18] and, in this respect, can 
be called quasi2D [20]. Kinematically the gas is two-dimensional, but the value of an effective 
2D coupling constant g for the interparticle interaction depends on the particle motion in the 
tightly confined direction. 

Let us first make a qualitative analysis of the interactions in such quasi2D gas. It can be 
viewed as a 3D gas which is uniform in two directions (x and y), and is confined to zero point 
oscillations by a harmonic potential muj^z 2 /2 in the third direction (z). The quasi2D regime 
requires the inequality 

tuu > n\g\ , (2.5) 

where n is the number of particles per unit area in the x — y plane and represents the 2D density 
of the quasi2D gas. Then the distribution of the density in the z direction is not influenced by 
the interactions and is given by 



n 3o(z) = — =exp -72 > ( 2 - 6 ) 




with lo = (h/mcoo) 1 ^ 2 being the harmonic oscillator length. The average interaction energy per 
particle is obtained by averaging the 3D interaction over the density profile in the z direction: 

T 93Dl^ri 2 3D (z)dz 
J_ oo n 3D {z)dz 

Eqs. H2.6fl and (|2.7|l lead to the following relation for the effective coupling constant: 

^2 



2V2H 2 a 



(2- 



171 Iq 

Accordingly, the criterion of the weakly interacting regime given by Eq. (|2.3|l takes the form 

/o»M • (2.9) 

One can easily check that under conditions (|2.5|l and (|2.9|l the gas satisfies the 3D criterion of 
weak interactions, ri3£i|a| 3 <C 1. 

The criterion (|2.9|) is independent of the gas density and in this respect is different from the 
criterion of the weakly interacting regime in the purely 2D case, following from Eqs. (|2.3|) and 
H2.4fl . The density enters the problem only through the condition of the quasi2D regime given 
by Eq. I|2.5|1 . However, the above analysis does not take into account the 2D character of the 
particle motion at large distances between them. In the next section we discuss the quasi2D 
scattering problem and derive a more exact criterion of weak interactions in the quasi2D case. 
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2.2 Quasi2D scattering problem 

Collisional properties of cold atoms strongly confined in one direction have been of great 
interest in the studies of spin-polarized atomic hydrogen. The interest was related to inelastic 
and elastic collisions in the 2D gas of hydrogen atoms adsorbed on liquid helium surface (see [21] 
for review). The creation of atomic gases in the quasi2D regime [16-19] opens new handles 
on studying 2D features of interparticle collisions, and provides spectacular evidence for the 
relation between the quasi2D scattering parameters and those in 3D. Theoretical studies of 
elastic and inelastic interactions in the quasi2D regime have been performed in [20,22,23], and 
it has been shown how a decrease of the tight confinement transforms this regime into the 3D 
one [22]. 

Here we consider elastic scattering of two atoms in the quasi2D regime [22]. So, the atoms 
are (tightly) confined to zero point oscillations in the axial (z) direction, and their motion in 
two other (x, y) directions is free at a large separation in the x — y plane. For a harmonic axial 
confinement, the motion of two atoms interacting with each other via the potential U(r) can 
be still separated into their relative and center-of-mass motion. The latter drops out of the 
scattering problem. The relative motion is governed by the potential U(r) and by the potential 
Vh{z) = itiujqZ 2 /4 originating from the axial confinement with frequency u>o- In the quasi2D 
regime one has the condition 

fiw„»£ , (2.10) 

where e = h 2 q 2 /m and q are the energy and wavevector of the motion in the x — y plane. 
Therefore, the atoms are in the ground state of the potential Vh(z) both in the incident and 
the scattered wave. The wavefunction of the relative motion satisfies the Schrodinger equation 

.^A + tf(r) + VffC0-^W)=^(r) • (2.11) 

We will consider the ultracold limit where the characteristic de Broglie wavelength of atoms 
greatly exceeds the radius of interatomic interaction R e . For the motion in the x — y plane the 
de Broglie wavelength is ~ 1/q, and we immediately obtain the inequality qR e <C 1. Under 
this condition the scattering is determined by the contribution of the s-wave for the motion in 
the x — y plane. In the axial direction, the atoms are tightly confined and the axial harmonic 
oscillator length Iq plays the role of their axial de Broglie wavelength. Therefore, we also 
require the condition 

h » Re , (2.12) 

which will allow us to consider only the s-wave for the three-dimensional relative motion of 
the atoms when they approach each other to short distances. The condition of the quasi2D 
regime (I2.10|) can be also written as qlo <C 1, and one clearly sees that Eqs. p.lOjl and (|2.12|l 
automatically lead to the inequality qR e <C 1. 

The scattering amplitude is defined through the asymptotic form of the wavefunction ip at 
an infinite separation p in the x — y plane, where it is represented as a superposition of an 
incident and scattered wave: 



V(r) « po(z) \e^P - f(q, ^ g^^'J , (2-13) 

where (fo(z) = (1/2-kI 2 ) 1 / 4 exp(— z 2 /4/q) is the eigenfunction of the ground state in the potential 
Vh{z), and <p is the scattering angle. The s-wave scattering is circularly symmetric in the x — y 
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plane, and the scattering amplitude is independent of <j>. Note that f(q,(j>) in Eg. (|2. 13(1 is 
different by a factor of —y/ETrq from the definition of the 2D scattering amplitude used in [24]. 

Relying on the condition l|2.12|) we will express the scattering amplitude through the 3D 
scattering length. At interparticle distances r 3> R e the relative motion in the x — y plane is 
free, and the motion along the z axis is governed only by the harmonic potential Vh(z). Then, 
the solution of Ea. (|2.11|l with U(r) =0 can be expressed through the Green function G £ (r, r') 
of this equation. Retaining only the s-wave for the motion in the x — y plane, we have 

f/>(r)=Mz)Mqp)-fGeir,0)/<j>o(fi) • (2-14) 



For p — > oo the Green function is G(r, 0) = <fio(z)4>o(0)y (i/Snqp) exp (iqp) and Ea. (|2.14|l gives 
the s-wave of Eq. (j2.13jl . 

The condition Iq 3> R e ensures that the relative motion of atoms in the region of interatomic 
interaction is not influenced by the axial (tight) confinement. Therefore, the wavefunction ip( r ) 
at distances R e <C r <C lo differs only by a normalization coefficient from the 3D wavefunction 
of free motion at zero energy, ip3D( r )- Writing this coefficient as r]ip o (0) and recalling that for 
r ^> R e one has ip3D(r) = (1 — a/r), we obtain 

ip(r) « Wo (0) (l - ~) ; i? e «r«; ■ (2.15) 

Ea. l2.15|) contains the 3D scattering length a and serves as a boundary condition for ?/>(r) 
(jTHjl at r -> 0. 

For r — > 0, a straightforward calculation of the Green function G(r,0) yields [20,22] 



G.(r,0)«-L- ' 



Attt 2(2tt) 3 / 2 Zo 



(2.16) 



where B w 0.915. With the Green function (|2.16|) . the wavefunction (|2.14|l at r — > should 
coincide with ip(r) l|2.15|l . This gives the coefficient rj and the scattering amplitude: 



/(e) 47r^ (0)ary ^ + (l/V2^)[ln(BfkJ /Tts) + m] ' 

One can see from Eq. I|2.17|l that the scattering amplitude is a universal function of the param- 
eters a/lo and e/hu>o- 

The 2D kinematics of the relative motion at interatomic distances p > lo manifests itself in 
the appearance of logarithmic dependence of the scattering amplitude on ejlhiOQ. The quasi2D 
scattering amplitude can be represented in the purely 2D form obtained, for example, in [24]: 

/(g) = hr(l/Jj +W 2 ' (2 - 18) 

In the purely 2D case a characteristic length d* depends on the shape of the interatomic 
potential U(r). For the considered quasi2D regime this length is expressed through Iq and the 
3D scattering length a: 

d q / 7r ; / / 7r ^o 



, e =V^ oexi VV2^ ' (2 - i9) 

where C is the Euler's constant. 

We thus conclude that the scattering problem in the quasi2D regime is equivalent to the 
scattering in an effective purely 2D potential which leads to the same value of d* . For a negative 
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a satisfying the condition \a\ <C lo, the effective potential is a shallow well which has a depth 
| Vo| and a radius Iq. This shallow well supports a weakly bound state with an exponentially 
small binding energy Eq, which leads to an exponentially large as follows from Ea. (|2.19[l . 
As a result, we have a resonance energy dependence of the scattering amplitude / at a fixed 
ratio a/lo, and a resonance behavior of / as a function of a/lo at a fixed e/hwo- 

Equation (|2.17|l shows that the resonance is achieved at a = a* (e) = —\Z2ttIq / In [BfvjjQ / ne). 
One can think of observing the resonance dependence of / on a/lo by measuring the rate of 
elastic collisions, which is proportional to |/| 2 . For example, one can keep temperature and too 
constant and vary a by using Feshbach resonances. This will be a striking difference from the 
3D case, where the rate constant of elastic collisions monotonically increases with a 2 . 

Stanford [25, 26] and ENS [27, 28] experiments with a thermal gas of Cs atoms tightly 
confined in one direction, observed a pronounced deviation of collisional rates from the 3D 
behavior. In these experiments the temperature was of the order of Hloq and, in this respect, 
they were in between the quasi2D and 3D regimes. 

The scattering amplitude / determines the coupling constant g for a fixed collision energy 
e. Away from the resonance, one may omit the imaginary part in Ea. (|2.17|l and / becomes 
real. Using Eg. (|2. 11(1 one can establish that in the ultracold limit the interaction energy for a 
pair of atoms is equal to g/£l, where the coupling constant is given by 

g = fd 3 r^ (z)U(r)^(r) , (2.20) 

and Q is the volume of the system. Accordingly, the total interaction energy is equal to (g/£l) 
multiplied by the number of pairs N 2 /2, as stated in the beginning of this Lecture. 

The main contribution to the integral in Ea. l|2.2U|l comes from interparticle distances r ~ 
R e . Therefore, one may put <poi z ) ~ Vo(0) and ip( r ) ~ Wo (0)^3-0 (?")■ Then, using the 
well-known result Jd 3 r ip3D(r)U(r) = Amf? a / m and Ea. l|2.17jl . we obtain 

g= *m = ^ i . (2 . 21) 

m m ^/27rl /a + ]xi(BhLJ /iTs) 

As one can see, the coupling constant for the interaction between particles depends on their 
relative energy. In a thermal gas, averaging the coupling constant over the energy distribution 
of particles simply leads to the replacement of e by T in Eq. I|2.21|) . 

In a Bose-condensed 2D gas, similarly to the 3D case (see, e.g. [29]), to zero order in 
perturbation theory the coupling constant is equal to the vertex of interparticle interaction in 
vacuum at zero momenta and frequency Q = 2/i. For low to > this vertex coincides with 
the amplitude of scattering at relative energy Co and, hence, is given by Ea. l|2.5|l with e = 2/i. 
This certainly requires the quasi2D condition Hujq 3> /i. For fi < 0, analytical continuation of 
Ea. H2.21[) to e < leads to the replacement e — > |e| = 2|/x|. 

We can now analyze the criterion of the weakly interacting regime for the ultracold quasi2D 
gas, assuming that it is dilute and the condition nR 2 < 1 is fulfilled. Since the scattering 
amplitude / and, hence, the coupling constant g can be written in the purely 2D form (see 
Ea. l|2.18|l ). the weakly interacting regime requires the inequality m\g\/2irh 2 <C 1 given above 
by Ea. H2.3fl . For a > 0, the quasi2D resonance is absent and this inequality is satisfied for 
any ratio a/lo. For a negative a, the system should be far away from the resonance and one 
should use Ea. (|2.21|l to make sure that the condition (|2.3|l is satisfied. This is demonstrated 
in Fig. 3 where we present the parameter mg /2irh 2 as a function of the ratio a/lo- The results 
are obtained from Ea. (|2.21JI for a fixed ratio fi-wo/e = 10 3 . In the vicinity of the resonance, 
where the criterion m\g\/2Trh 2 <g; 1 is not satisfied, they are shown by the dashed curves. 
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Fig. 3. The parameter mg/2nh 2 versus a/lo, obtained from Ea. l|2,21|l (see text). 



2.3 True BEC at T = 

We start the discussion of BEC regimes in (quasi)2D weakly interacting gases with the case of 
zero temperature. From this point on we consider only repulsive interaction between particles 
(g > 0). The behavior of attractively interacting Bose gases and, in particular, the problem of 
collapses are beyond the scope of these lectures. 

In the second quantization the Hamiltonian of the system reads: 

H = jdv¥ (-^A + V(r) + , (2.22) 

where ^(r, t) is the Heisenberg field operator of the atoms, V(r) is an external (trapping) 
potential, and all field operators in the integrand are taken at the same position r and for the 
same time t. The Heisenberg equation of motion for the field operator takes the form: 

ih— = [*,#]= (-— A + Vtf+gtf*)* . (2.23) 

We now assume a priori that there is a true condensate at T = 0, and the condensate density 
no is much larger than the density of non-condensed particles n'. Accordingly, we represent 
the field operator as a sum of the non-condensed part iff' and the condensate wavefunction 

= \fno which is a c-number: 

* = + * ■ (2.24) 

At equilibrium, the time dependence of the condensate wavefunction is reduced to <x 
exp (—ifit/H), where fi is the chemical potential. Then, taking average of both sides of Ea. H2.23|) 
and omitting the contribution of non-condensed particles, we obtain the Gross-Pitaevskii equa- 
tion for the condensate wavefunction [30,31]: 

|^A + U(r)+ 5 |*o| 2 -M)*o-0 . (2.25) 

The function is normalized by the condition 

ldr\^ \ 2 = N , (2.26) 



WILL BE SET BY THE PUBLISHER 



Prl-13 



which gives a relation between the number of particles N and chemical potential /x. 

Equations l|2.23[l and (|2.25[1 immediately lead to the equation of motion for the non- 
condensed part of the field operator: 

ih—=(- — A + V(r) + 2g\* \ 2 )*' + gt> 2 Q ^ . (2.27) 



dt \ 2m 

We then use the Bogoliubov transformation (see [32]) generalized for the spatially non- uniform 
case [33] and express through the eigenmodes of elementary excitations: 

ty' = cxp(—ifit/h) ^ \u v a v exp(— ie v t/K) — v*al exp(ie„t/h)] . (2.28) 

V 

Here a U) aj, are annihilation and creation operators of the excitations, and £„ are their eigenen- 
ergies. Excitations are characterized by a set of quantum numbers v, and their Bogoliubov 
{u„,t>„} functions are normalized by the condition 

dr (K(r)| 2 - K(r)| 2 ) =1 . (2.29) 

Commuting both sides of Eg. 1)2.27(1 with the operators a v and a J,, we arrive at the Bogoliubov- 
de Gennes equations for the excitation energies and wavefunctions. Assuming that aside from 
the factor exp(—ipit/h), the condensate wavefunction is real, these equations read: 

-^A + V(r)+g\y \ 2 -^f+ = e v f- (2.30) 

-^-A + V(r) + 3g\y a \ 2 - f ?j f- = e u f+ , (2.31) 

where the functions /jf =u v ±v v . 

In the spatially uniform case Ea. (|2.25|l is reduced to the expression for the chemical potential 
through the condensate density: [i = nog. Excitations are characterized by their wave vector 
k, and Eqs. I|2.30|l . (|2.31() give the eigenfunctions 



j / \ ±1/2 

/± = u k ± v k = - ( ^- ) exp (ik.r) (2.32) 



and lead to the Bogoliubov energy spectrum 



£ fc = V E t + 2 ^Ek , (2.33) 

where Ek — h 2 k 2 /2m is the energy of a free particle. For small momenta, excitations are 
phonons characterized by a linear dispersion law Ek = hc s k, where c s = U fi/m is the velocity 
of sound. High-momentum excitations are single particles with ek = Ek + \i. Characteristic 
momenta at which the linear spectrum transforms into the quadratic one are of the order of 
l// c , where l c — h/ y/mpt is called correlation or healing length. The corresponding excitation 
energies are ~ /i. Note that the criterion l|2.3|) of the weakly interacting regime is equivalent 
to the condition l c ^S> f. In other words, there are many particles under the healing length. 

We now calculate the density of non-condensed particles n' and the single-particle correla- 
tion function. Using Eqs. (|2.28|) and (|2.32|l we obtain 

n' = (*Hr)*\r)) = J^vi = ^n . (2.34) 




One clearly sees that the ratio n! /no is small. It is proportional to the small parameter mg/2irh 2 
required for the weakly interacting regime in the 2D Bose gas. 

The single-particle correlation function is g\ (r', r") = (^'(r')^f(r")) and it depends only on 
the relative coordinate r = |r' — r"|. On the basis of Eqs. (|2.24|) . 12.28|) . and H2.32[) we obtain 



9i (r) =n + 



nd 2 k 



v 2 exp(ik.r) = no 



1 



mg I r 

2^? 1 r 



Ki 



(2.35) 



where I\ and K\ are the growing and decaying modified Bcssel functions, respectively. The 
result of Ea. (|2.35l) is displayed in Fig. 4. For r -C l c the single-particle correlation function is 
close to the total density n. At distances r 3> l c , Ea. H2.35|l gives 



9i(r) = n 1 + 



mg l c 
Anh 2 r 



r > lc 



(2.36) 



For r — > oo the second term in the brackets vanishes and g\ is tending to no- Thus, there is a 
long-range order in the system. This justifies our initial assumption of the presence of a true 
BEC in the 2D weakly interacting gas at T = 0. Note that g\ drops from n to no at distances 
of the order of the healing length l c . 

In order to go beyond the Bogoliubov approach and find, for example, corrections to the 
excitation energies, which are small as mg/2irh 2 , one may proceed along the lines of the Bcliaev 
theory [34] developed for 3D Bose-condensed gases at T = 0. 



2.4 Quasicondensate at finite temperatures 

As we already mentioned, the finite-temperature uniform 2D Bose gas is characterized by the 
absence of a true Bose-Einstein condensate and long-range order. To gain insight into the 
nature of this phenomenon we first try to calculate the density of non-condensed particles, 
assuming that there is a true BEC and one can use the Bogoliubov approach described in the 
previous subsection. Then, at finite temperatures, from Eqs. (|2.28(l and (|2.32() we obtain 



n'= (tf't(r)tf'(r)) =n' v + n' T 



(2.37) 
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The vacuum (temperature- independent) contribution n' v is given by Eg. (|2. 34(1 . and the thermal 
contribution is determined by the relation 

n<fk , 9 9 , f d 2 k ( E k + n\ 

(v 2 + ul)N k = / —- ^—^ N k , (2.38) 



(2tt) 2 KJ * J (2tt)2 V e k 

where N k — [exTp(e k /T) — l]^ 1 are equilibrium occupation numbers for the excitations. For 
k — » the occupation numbers are N k ss T/e k , and the excitation energy is e k cx k. Therefore, 
the integrand in the rhs of Ea. H2.38[l behaves as d 2 k/k 2 and the integral is divergent at small 
momenta, which rules out the assumption of the presence of a true condensate. The origin 
of this infrared divergence is related to long-wave fluctuations of the phase of the condensate 
wavefunction, which can be understood turning to the density-phase representation for the 
field operators: 

VP = cxp(i(j>)Vfi, # f = V^exp(-i^) . (2.39) 
Here the density and phase operators are real and satisfy the commutation relation 

[n(r), 4>(r')} = iS(r - r') . (2.40) 

We now present a general approach which is based on the Hamiltonian l|2.22[) and includes 
the presence of an external trapping potential V(r). We assume a priori and justify later that 
fluctuations of the density are small. Substituting Eqs. I|2.39|l into Eq. I|2.23|l and separating 
real and imaginary parts, we get the coupled continuity and Euler hydrodynamic equations for 
the density and velocity v = (7i/m)V<^: 

-h-^— = — (A0V^ + 2V0V\/^) , (2.41) 
at 2m 

dd> h 2 h 2 

-h^-Vft = — (V0) 2 \/^- — &Vfi + V(Y)Vfi + gh 312 . (2.42) 
at 2m 2m 

For small density fluctuations, Ea. H2.41|) shows that fluctuations of the phase gradient are 
also small. Writing the density operator as n(r) = n(r)+5n(r) and shifting the phase by —/j,t/H, 
we then linearize Eqs. I|2. 4112.42(1 with respect to Sh, around the stationary solution h = n, 

= 0. The zero order terms give the Gross-Pitaevskii equation for the mean density n: 

' n +V(r)+gn = n , (2.43) 



2m y 

and the first order terms provide equations for the density and phase fluctuations: 

hd[6n/y/n)/dt = (-h 2 V 2 /2m + V(r) + gn - j u)(2 N /n0) , (2.44) 
-~hd(2y/^^)/dt = (-h 2 V 2 /2m + V{r) +3gn- (i)(5h/^) . (2.45) 

Solutions of Eqs. (|2.44I2.45|I for 5n, V(j> are obtained in terms of elementary excitations: 

5n(r) = n(Y) 1/2 ^f-(v)e- l£ » t/h a v + \i.c. , (2.46) 

0(r) = [4n(r)]- 1 /^_ i/ +( r ) e --,tAa i , + h.c. , (2.47) 

V 

where the eigenfunctions f„ obey the Bogoliubov-de Gennes equations l(2.3t)|l and (|2.31|l , with 
l^o 1 2 replaced by n. In the previous subsection these functions were introduced as eigen- 
functions of elementary excitations of a true Bose-Einstein condensate. We thus see that the 
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assumption of small density fluctuations is sufficient for having the Bogoliubov wavefunctions 
and spectrum of the excitations, irrespective of the presence or absence of a true conden- 
sate. Note that this statement holds for both uniform and trapped Bose gases in the weakly 
interacting regime in any dimension. 

In the uniform case the chemical potential is y, = gn, and the excitation wavefunctions and 
spectrum are given by Eqs. I|2.32[l and l|2.33[) . As discussed above, the characteristic excitation 
energies at which the character of the spectrum transforms from the phonon into the single- 
particle one are of the order of \i. We, therefore, separate the energy space into two regions: 
low-energy (phonon) part with Ek < n, and high-energy (free-particle) part where Ek > [i- The 
operators of the density and phase fluctuations are then represented as: 

Sf\ = Sfip + Shf , (2.48) 
$ = + , (2.49) 

where the indices p and / stand for the phonon (e < /i) and free-particle (e > /i) parts, 
respectively. 

Fluctuations originating from the high-energy part are small. This is seen from the cal- 
culation of the density of this part of the gas, < $\$f >, where the operator 4ff ac- 
counts for both the density and phase fluctuations and is given by the Bogoliubov trans- 
formation H2.28J) in which the summation is performed only over excitations with energies 
Ek > /U. At energies significantly larger than \i the function Vk = (f£ — f^)/2 — > 0, and 
u k = (fk + fk )/^ ~ * n -1 ' 2 exp(ik • r). The excitation energy is ~ E(k) + fj, and, hence, 
the operator ^ / describes an ideal thermal gas of Bose particles with chemical potential equal 
to —fj,. The density of this gas is exponentially small at T < [i. In the 2D case, for higher 
temperatures it is equal to 



and is much smaller than n assuming that T is well below the temperature of quantum degen- 
eracy Td = 2irh 2 n/m. 

The low-energy fluctuations of the density at T <C Td are also small. This follows from 
the calculation of the density-density correlation function (Sn s (r)8h s (0)). For the 2D gas, a 
straightforward calculation using Ea. l|2.46|l yields 

= ^ £ lU (2Nk + 1} C ° S(k - r) < ^ Wd^ } 1 • (2 ' 50) 

Similarly, one finds that the density-phase correlation functions ^"p( r ) <ft p (0)) and (<ftp(r) a ) 
are small for T <C T d . 

To zero order in perturbation theory, omitting small high-energy fluctuations and small 
low-energy fluctuations of the density, the single-particle correlation function is found by using 
the field operator in the form: 

* = V^exp(i0 p ) . (2-51) 

Relying on the Taylor expansion of the exponent one proves directly that (exp{i0 p }) = 
exp(— (<fip}/2). Then, for the single-particle correlation function we obtain: 



9l (r) = ($t(r)#(0)) = ncxp f-i<(4(0) - ^ p (r)) 2 ) 



T < T d . (2.52) 
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In the 2D case, vacuum low-energy (long-wave) fluctuations of the phase are small as 
mg/2-Kh 2 . However, thermal phase fluctuations are large for r — > oo. On the basis of Eg . 1)2 .47(1 . 
we find the following asymptotic expression for the mean square phase fluctuations: 

p (r)-0 p (O)] 2 ^«^ln^^) ; r » A T . (2.53) 

where At is equal to the healing length l c at T ^> fx, and to the thermal de Broglie wavelength 
of phonons hc s /T for T <C /i. Accordingly, the correlation function <?i(r) undergoes a slow 
power law decay at large distances: 

/A \ T/Td 

9i( r )= n [-f) ; ^»At ■ (2.54) 

This is drastically different from the situation in 3D Bose-condensed gases, where fluctuations 
are small at any distance and g\ is tending to the condensate density as r — ► oo. We thus see 
that just long-wave thermal fluctuations of the phase destroy the long-range order and true 
BEC in finite-temperature 2D Bose gases. 

The low-temperature behavior of the single-particle correlation function, showing its power 
law decay at large r, was first obtained by Kane and Kadanoff [9], The idea of dividing the 
system into slow (low-energy) and fast (high-energy) parts belongs to Popov who developed 
a perturbation theory for Bose systems on these grounds [1]. In the described hydrodynamic 
approach, for obtaining perturbative corrections to Eq. I|2.52|) one should expand the continuity 
and Euler equations l|2.4I|l and l|2.42|l up to second order in the density and phase fluctuations, 
which provides a correction to the stationary solution. One should then include the low-energy 
density fluctuations and the high-energy fluctuations in the expression for the field operator. 
An expression for g\(r) in the uniform 2D case, which contains these corrections and is obtained 
on the basis of the Popov theory, is given in [35]. 

The distance at which the mean square phase fluctuations become of the order of unity and 
the single-particle correlation function significantly decreases, is called the phase coherence 
length 1$. From Eq. (|2.53|l we obtain 

^ = A T expQj^ (2.55) 

and clearly see that at T <C the phase coherence length greatly exceeds the healing length 
l c . This means that the system can be divided into blocks of size L which is chosen such that 
l c <C L <C V Then, using Eqs. (|2.52|) and (|2.53|) we can make sure that correlation properties 
inside each block are the same as in a genuine Bose-condensed gas. We thus conclude that there 
is a true condensate in each block. However, the phases of different blocks are not correlated 
with each other. Therefore, this system got the term quasicondensate, or condensate with 
fluctuating phase [13]. For the existence of the quasicondensate, it is crucial that the density 
fluctuations are small on any distance scale. 



2.5 True and quasicondensates in 2D traps 

Bose-Einstein condensation in trapped 2D gases is qualitatively different from that in the 
uniform case. As was discussed in Lecture 1, for an ideal 2D gas the change of the density of 
states due to a harmonic confining potential, leads to a macroscopic occupation of the ground 
state of the trap (ordinary BEC) at temperatures T < T c paJV 1 ' 2 ^, where TV is the number 
of particles, and uj the trap frequency [7]. Thus, there is a question of whether an interacting 
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trapped 2D gas supports the ordinary BEC or the Kosterlitz-Thouless type of a cross-over to 
the BEC regime. Related studies are now underway, and for large 2D samples one expects the 
Kosterlitz-Thouless cross-over. However, irrespective of the type of the BEC cross-over, the 
critical temperature will be always comparable with T c of an ideal gas. On approach to T c 
from above, the gas density is n c ~ iV/i? 2 ^, where Rt c ~ y/T c /muj 2 is the thermal size of the 
cloud, and hence the Kosterlitz-Thouless temperature is ~ h 2 n c 

An important feature of trapped gases is that the confining potential introduces a finite 
size of the sample, which sets a lower bound for the momentum of elementary excitations 
and reduces the phase fluctuations. For this reason, at finite temperatures well below T c the 
phase fluctuations are small, and the equilibrium state is a true condensate. At intermediate 
temperatures T < T c the phase fluctuates on a distance scale smaller than the size of the gas 
sample, and one has a quasicondensate. Qualitatively, the character of the BEC state can be 
identified by comparing the size of the sample with the uniform-gas phase coherence length 
at the maximum density of the trapped gas. In this Lecture we present a more detailed 
analysis of the phase fluctuations and BEC character in weakly interacting trapped (quasi) 2D 
gases. For simplicity we assume that the coupling constant g is independent of the gas density. 
In the quasi2D case this requires the condition Iq ^> \a\ which leads to g given by Ea. l|2.8|l . 
This condition is satisfied in recent experiments [16-19], where the quasi2D regime was reached 
for atomic Bose condensates. 

For finding the phase (and density) fluctuations in a trapped gas, one should know the 
density profile and the spectrum of elementary excitations. If the number of particles is very 
large and the chemical potential fi greatly exceeds the level spacing in the trap, the kinetic 
energy term in Eq. (|2.43|1 is much smaller than the nonlinear term and can be neglected. This 
approach is called the Thomas-Fermi (TF) approximation, and in a harmonic trap the density 
profile takes the well-known parabolic shape: 




(2.56) 



with u>i and r» being the trap frequency and coordinate in the i-th direction. The dependence 
of the excitation spectrum on the trapping geometry has been extensively studied for 3D TF 
condensates (see [39] for review) . The spectrum and wavefunctions of low-energy excitations 
(e„ <C /i) can be found analytically [36-38]. For the ID and 2D weakly interacting trapped 
Bose gases this has been done in [40,41]. 

We analyze the character of BEC in a harmonically trapped symmetric 2D gas with repulsive 
interparticle interaction, relying on the calculation of the single-particle correlation function 
[20]. As well as in the uniform case, fluctuations of the density are provided by excitations 
with wavelengths of the order of l/l c and in a similar way one proves that they are small 
at temperatures T -C T^. Therefore, the single-particle correlation function can be found by 
using the field operator from Ea. l2.51|) . with coordinate-dependent mean density n(r) and 
phase operator 4> p {r). Then the correlation function is given by Ea. (|2.52Jl generalized to the 
trapped case: 

(*t (r )^( )) = VMr)n (0)exp ~ M*)?)) , (2-57) 

where r = at the trap center. 

In the Thomas-Fermi regime, the radius of the gas sample is Rtf = (2/i/mo; 2 ) 1 / 2 , and 
integration of Eq. (|2.56|) over the spatial volume of the gas gives a relation between the chemical 
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potential and the number of particles: 



I Nmg 




(2.58) 



/! = n 0m g = 



where n m = n(0) is the maximum density. It is worth noting that the ratio of the chemical 
potential to the BEC transition temperature T c of an ideal gas is fi/T c ~ (rag/irk 2 ) 1 ^ 2 <C 1. 

For calculating the mean square fluctuations of the phase, we use the (discrete) spectrum 
and wavefunctions of excitations with energies e v -C /i, obtained relying on the method devel- 
oped for Thomas- Fermi trapped condensates [36-38,40,41]. At temperatures larger than the 
chemical potential, for distances r greatly exceeding the healing length l c at the trap center, 
we obtain a result similar to that in the uniform case: 



The character of the Bose-condensed state is determined by the phase fluctuations at dis- 
tances r ~ Rtf- If they are small, one has a true condensate, and for a large value of these 
fluctuations the state is a quasicondensate. With logarithmic accuracy, from Ea. (|2.59|l we find 



In the case of a quasicondensate the phase coherence length 1$ is given by an expres- 
sion similar to that for uniform quasicondensates. For T > fi, equation (|2.59(l yields 1$ « 
l c ex-p(irh 2 n m /mT). As 1$ greatly exceeds the healing length, the quasicondensate has the 
same Thomas-Fermi density profile as the true condensate. Correlation properties at dis- 
tances smaller than 1$ and, in particular, local density correlators are also the same. Hence, 
one expects the same reduction of inelastic decay rates as in 3D condensates [13]. However, 
the phase coherence properties of a quasicondensate are drastically different. We will have a 
detailed discussion of this subject in the next Lecture for the case of ID quasicondensates. 

From Eq. (|2.60[1 we see that in large gas samples one has more possibilities for obtaining the 
quasicondensate regime. In quasi2D trapped alkali gases one can expect the number of trapped 
atoms N ranging from 10 4 to 10 6 and a value ~ 10~ 2 or larger for the small parameter mg/2irh 2 . 
In the MIT sodium experiment [16] on achieving the quasi2D regime for a trapped condensate, 
this small parameter was about 1CP 2 and the number of atoms was ~ 10 5 . Therefore, even at 
temperatures somewhat higher than fi the gas was in the regime of true BEC. Note, however, 
that the coupling constant g in the MIT experiment [16] was rather small as the scattering 
length for sodium is only 28 A. One can think of achieving the quasicondensate regime for the 
same trapping parameters and N as at MIT, by using Feshbach resonances and tuning the 
coupling constant to much larger values. 

Presently, the largest obtained 2D Bose-condensed gas is the one of spin-polarized atomic 
hydrogen on liquid helium surface [15]. In this system the temperature is about 100 mK, the 
density exceeds 10 13 cm~ 3 , and the number of particles is ~ 10 s . Estimates show that the 
state should be a quasicondensate rather than a true condensate. However, this should still be 
proven experimentally. 

Problem: Prove that to zero order in perturbation theory the coupling constant for the 2D 
Bose-condensed gas is g — h 2 f/m, where / is the scattering amplitude (as defined in Ea. (j2.13|i l 
at energy of the relative motion E = 2/i. (Yu.E. Lozovik, 1971). 



(A(o)-^ S (%))'«pJin(:) 



(2.59) 




(2.60) 
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3 LECTURE 3. TRUE AND QUASICONDEISISATES IN ID TRAPPED GASES 

One-dimensional Bose systems at low temperature show a remarkable physics not encountered 
in 2D and 3D. In particular, the ID Bose gas with repulsive interparticle interaction becomes 
more non-ideal with decreasing ID density n [42,43]. The regime of a weakly interacting gas 
requires the correlation length l c = h/y/mng to be much larger than the mean interparticle 
separation 1/n. For small n or large interaction, where this condition is violated, the gas 
acquires fermionic properties as the wavefunction strongly decreases at short interparticle dis- 
tances [42,43]. In this case it is called a gas of impenetrable bosons or Tonks- Girardeau gas 
(cf. [44]). 

Spatially uniform ID Bose gases with repulsive interparticle interaction have been exten- 
sively studied in the last decades. For the delta-functional interaction, Lieb and Liniger [43] 
have calculated the ground state energy and the excitation spectrum which at low momenta 
turns out to be phonon-likc. Generalizing the Lieb-Liniger approach, Yang and Yang [45] 
have proved the analyticity of thermodynamical functions at any finite temperature T, which 
indicates the absence of a phase transition. 

However, at sufficiently low T the correlation properties of a ID Bose gas are qualitatively 
different from classical high-T properties. In the regime of a weakly interacting gas the density 
fluctuations are suppressed [9] , whereas at finite T the long- wave fluctuations of the phase lead 
to an exponential decay of the single-particle correlation function at large distances [9,46]. A 
similar picture, with a power-law decay of the density matrix, was found at T — [47,48]. 
Therefore, the Bose-Einstein condensate is absent at any T, including T — 0. These findings 
are consistent with the Bogoliubov fc~ 2 theorem at finite T, and with a related treatment at 
T = [49]. 

Earlier studies of ID Bose systems (see [1] for review) allow us to conclude that in uniform 
ID gases the decrease of temperature leads to a continuous transformation of correlation prop- 
erties from ideal-gas classical to interaction/statistics dominated. In the weakly interacting 
regime at low T, where the density fluctuations are suppressed and the phase fluctuates on a 
distance scale much larger than the correlation length l c , one can speak of a phase-fluctuating 
Bose-condensed state (quasicondensate) . A ID classical field model for calculating correlation 
functions in the conditions, where both the density and phase fluctuations are important, was 
developed in [50] and for Bose gases in [51]. The use of the Bogoliubov approach for describing 
ID uniform (quasi)condensates is recently discussed in [52] and is presented in the lectures of 
Yvan Cast in. 

In order to avoid misunderstanding, we should make a remark on the presence of a qua- 
sicondensate in the uniform weakly interacting ID Bose gas. It is sometimes stated that in 
this case the quasicondensate exists only at T = 0, and its presence is related to the power 
law decay of the single-particle correlation function at large distances (see, for example, [53]). 
However, at finite temperatures where this correlation function decays exponentially, the phase 
coherence length can greatly exceed the healing length l c . This is the case for small density 
fluctuations, and then the physical picture is similar to that in uniform finite-temperature 2D 
gases or zero-temperature ID gases. The system can be divided into blocks of size L satisfying 
the inequalities I c « L « 1^, and inside each block one has a true condensate with a phase 
that is not correlated with phases in other blocks. In this sense, the term quasicondensate can 
also be used in weakly interacting uniform ID Bose gases at finite temperatures. 

The ID regime for trapped atomic condensates has been achieved in several experiments 
by tightly confining the motion of particles in two directions [16,54,55]. In this Lecture we 
will discuss the various quantum degenerate states that are present in finite-temperature ID 
trapped Bose gases and focus attention on the case of a weakly interacting gas. 
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3.1 Weakly interacting regime in ID 

We will discuss ID Bose gases with repulsive short-range interaction between particles (coupling 
constant g > 0). Present realizations of the ID regime imply particles in a cylindrical trap, 
which are tightly confined in the radial direction, with the confinement frequency ujq greatly 
exceeding the mean-field interaction. Then, at sufficiently low T the radial motion of particles 
is essentially "frozen" and is governed by the ground-state wavefunction of the radial harmonic 
oscillator. If the radial extension of the wavefunction, Iq — (H/muJo) 1 / 2 , is much larger than 
the radius of the interatomic potential R e , the interaction between particles acquires a 3D 
character and is characterized by the 3D scattering length a. For this case, the coupling 
constant has been found by Olshanii [56] , and for 1 > 8 it can also be obtained by averaging 
the 3D interaction over the radial (harmonic oscillator) density profile: 

2h 2 a . . 

9= 72 5 l o > a . (3.1) 

m Iq 

Statistical properties of such quasilD samples are the same as those of a purely ID system 
with the coupling constant given by Ea. (|3.1|l . 

As well as in the 2D case, we obtain the criterion of the weakly interacting regime at T = 
by comparing the interaction energy per particle, I = ng, with a characteristic kinetic energy 
of particles at a mean separation f between them. In the ID case we have f ~ 1/n, and this 
kinetic energy is K sa h 2 n 2 /m. In the weakly interacting regime, where the wavefunction of 
particles is not influenced by the interaction at interparticle distances of the order of f , one 
should have / <C K . This leads to the criterion of the weakly interacting regime 

7=^«1 • (3-2) 

From Eq. l|3.2[l one really sees that in contrast to 2D and 3D gases, the ID Bose gas becomes 
more interacting with decreasing density. In the purely ID case this statement is valid as long 
as the mean interparticle separation 1 jn remains much larger than the radius of interparticle 
interaction R e . For quasilD weakly interacting gases realized in present experiments, the 
density should still be such that the healing length l c = h/^/mng ^> Iq, otherwise the gas 
leaves the ID regime. Assuming a <C Iq, this gives the condition iiflCl. 

The parameter 7 can also be interpreted as a ratio of the mean interparticle separation to the 
characteristic interaction length for two particles, r g = h 2 /rag. This quantity determines the 
distance scale on which the repulsion between the particles reduces their relative wavefunction. 
Under the condition r g 1/n, which is equivalent to Ea. H3.2(l . this reduction is practically 
absent and, hence, the gas is weakly interacting. In contrast, for r g <C 1/n or, equivalently, 
7> 1, the relative wavefunction is strongly reduced at distances smaller than 1/n. The gas 
then enters the strongly interacting Tonks- Girardeau regime and acquires fermionic properties. 

For particles trapped in a harmonic (axial) potential V(z) = mw 2 z 2 /2, one can introduce 
a complementary dimensionless quantity 

mgl . . 

a = —p- , (3.3) 



where I = y/hfmu) is the amplitude of axial zero point oscillations. The parameter a provides 
a relation between the interaction strength g and the trap frequency u. Actually, it represents 
the ratio of / to the interaction length r g . For a 1, the strongly interacting regime is not 
present at all as the relative motion of two particles on approach to each other is governed by 
their harmonic confinement rather than the interparticle interaction. 
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3.2 True BEC and diagram of states at T = 

The nature of quantum degenerate states in the trapped ID gas is strongly influenced by the 
interparticle interaction and by the presence of the trapping potential. The latter introduces 
a finite size of the sample and provides a low-momentum cut-off for the phase and density 
fluctuations. This reduces the fluctuations compared to the uniform case. 

We first discuss the case of T = 0. For simplicity, we will consider a uniform gas with a 
finite size L. We assume a priori that there is a true condensate with a density no greatly 
exceeding the density of non-condensed particles n' . Then, for a large value of L we have the 
Bogoliubov excitation spectrum given by Eq. I|2.33|> and the excitation wavefunctions following 
from Eq. l|2.32|l . The density of non-condensed particles is given by an equation similar to 
Eq. I|2.34|) • The difference is that now the lower limit of integration over the excitation momenta 
is non-zero. It is related to a finite size of the system and is equal to ir/L. We thus have 

nL = r*L± , ( 3.4) 

and a straightforward calculation yields 

no 7r \ ewl c J 

For a realistic size L of a trapped ID Bose gas, the logarithmic factor in Ea. H3.5l) does not 
exceed 10. Hence, for 7 <C 1 we have n' <C no. This proves that fluctuations are small and 
justifies our initial assumption on the presence of a true BEC in a weakly interacting trapped 
Bose gas at T = 0. An analysis of vacuum fluctuations for ID Bose gases in harmonic traps is 
given in [41]. 

Thus, in a harmonically trapped ID Bose gas at T = we have a true condensate with the 
density profile determined by the Gross-Pitaevskii equation l)2.25[l . In the Thomas- Fermi (TF) 
regime, where the chemical potential is /i» tko, the kinetic energy term can be omitted and 
we have a parabolic density distribution typical for condensates in harmonic traps: 

n (z) = n 0m ( 1 - ) (3.6) 

V H TF/ 

for —Rtf < z < Rtf, and is equal to zero otherwise. The (half) size of the condensate is 
Rtf = (2/i/niw 2 ) 1 / 2 , and the maximum density is no(0) = no m = fi/g. Integrating Ea. (|3.t)|) 
over dz we obtain a relation between the chemical potential and the number of particles N: 

where the parameter a is given by Ea. (|3.3|) . 

For a > 1 we are always in the TF regime. In this case, Ea. (|3.2f> requires a sufficiently 
large number of particles: 

N > N* = a 2 , (3.8) 

which reflects the fact that the weakly interacting regime requires sufficiently large densities. 
For a«l the criterion (|3.2() of a weakly interacting gas is satisfied at any N. The condensate 
is in the TF regime if N 3> a -1 . In the opposite limit the mean-field interaction is much 
smaller than the level spacing in the trap Hco. Hence, one has a macroscopic occupation of the 
ground state of the trap, i.e. there is an ideal gas condensate with a Gaussian density profile. 
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Fig. 5. Diagram of states for a trapped ID gas at T = 0. 



If a 3> 1 and N <C N*, then the trapped ID gas is in the strongly interacting Tonks- 
Girardeau regime. The one-to-one mapping of this system to a gas of free fermions [42] ensures 
the fermionic spectrum and density profile. The chemical potential is equal to Nfku, and the 
density distribution is n(z) — {s/2N — (z/R) 2 , where the size of the cloud is R — \/2Nl. 
Note that this profile is different from both the profile of the zero-temperature condensate and 
the classical distribution of particles. Correlation properties of strongly interacting ID Bose 
gases will be discussed in Lecture 4. 

In Fig. 5, we present the diagram of states for the zero-temperature trapped ID Bose gas in 
terms of N and a. This diagram clearly shows the presence of three states as discussed above: 
Thomas-Fermi BEC, condensate with a Gaussian density profile, and the Tonks-Girardeau gas. 

3.3 Regimes of quantum degeneracy at finite temperature 

The characteristic temperature of quantum degeneracy for a trapped ID Bose gas is Td ~ Nfuo. 
In the weakly interacting regime, the character of finite-temperature Bose-condensed states at 
T <C Td is determined by thermal fluctuations. Before describing these states, we briefly 
analyze a cross-over to the BEC regime in ID trapped gases. 

As we discussed in Lecture 1, for an ideal gas one has a sharp cross-over to BEC. Ket- 
terle and van Druten [8] found that the decrease of T to below T c = Nhui/ ln(2iV) strongly 
increases the population of the ground state, which rapidly becomes macroscopic. This sharp 
cross-over originates from the discrete structure of the trap levels. However, the presence 
of the interparticle interaction changes the picture drastically. One can distinguish between 
the (lowest) trap levels only if the interaction between particles occupying a particular level 
is much smaller than the level spacing. Otherwise the interparticle interaction smears out 
the discrete structure of the levels. For T close to T c the occupation of the ground state is 
~ T c /huj sw N/ ln(2iV) [8] and, hence, the mean-field interaction between the particles in this 
state (per particle) is Ng/l\n(2N). The sharp BEC cross-over requires this quantity to be 
much smaller than hu, and we arrive at the condition iV/ln(2iV) <C a -1 . For a realistic num- 
ber of particles (N ~ 10 3 — 10 4 ) this is practically equivalent to the condition N <C or 1 at 
which one has the ideal gas Gaussian condensate. 

As we see, the sharp BEC cross-over requires small a. For possible realizations of ID gases, 
using the coupling constant (|3.1J) for Iq 3> a, we obtain a — 2al/l 2 . Then, even for the ratio 
l/lo ~ 10 and moderate radial confinement with Iq ~ 1 /im, we have a ~ 0.1 for Rb atoms 
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(a ~ 50 A) . Clearly, for a reasonably large number of particles the (sharp) cross-over condition 
N <C or 1 can only be fulfilled at extremely small interparticle interaction. One can think of 
reducing a to below lA and achieving a < 10~ 3 by using Feshbach resonances. In this case 
one can expect the sharp BEC cross-over already for N ~ 10 3 . 

Otherwise, similarly to the uniform ID case, the decrease of temperature to below con- 
tinuously transforms a classical ID gas to the regime of quantum degeneracy. If the condition 
l|3.2[l is satisfied, one has a weakly interacting gas which at T = becomes a true TF con- 
densate. To understand the nature of a Bose-condensed state at finite T <C T<j, we analyze 
the behavior of the single-particle correlation function by calculating the fluctuations of the 
density and phase [57]. 

We assume a priori that the weakly interacting trapped ID Bose gas is characterized by 
small density fluctuations and justify this assumption later. Then the operators of the density 
and phase fluctuations are given by Eqs. H2.48J1 and l|2.49(l . and the spectrum and wavefunctions 
of elementary excitations follow from Eqs. H2.30[l and (|2.31|) in which the density distribution 
is the one of the zero-temperature condensate. It has already been shown in section 3.2 that 
vacuum fluctuations of the density and phase are small. Therefore, we now discuss only thermal 
fluctuations. Using Eq. l|3.7[) we find that fi/Td ~ (a 2 /A) 1 / 3 -C 1. Accordingly, our analysis 
will include both cases, T -C /i and T ^>> fi. 

Fluctuations coming from the high-energy part of the excitations (e v > fi) are small. 
Similarly to the uniform case discussed in Lecture 2, this part can be viewed as an ideal 
thermal gas of particles, with chemical potential equal to — /i. One then finds that the number 
of particles in the high-energy part becomes exponentially small for T <C fi, and for T ^> /j, 
it is ~ (T/T^N <C A. We therefore confine ourselves to fluctuations coming from low-energy 
excitations (e„ < //). 

The solution of Eqs. I|2.30|l and (|2.31|l for these excitations gives the spectrum e,j — 
hw\j + 1) [40,41] and wavefunctions 

ft = 



Rtf 



1/2 



Pj(x) , (3.9) 



where j is a positive integer, Pj are Legendre polynomials, and x — z/Rtf- Using Eqs. (|2.48|l 
and (|3.9|l . for the mean square (thermal) fluctuations of the density we have 

((5h(z) - 5h{z')Y) T ^ {5nl z ,) T = £in0m{ i + 1,2 \ p^) - P^')) 2 N 3 , (3.10) 

3=1 

with Nj — [exp(e J /T) — 1] being the occupation numbers for the excitations. Assuming T ^> Hcu, 
the main contribution to the density fluctuations comes from quasiclassical excitations (j 1). 
These fluctuations are largest on a distance scale \z — z'\ greatly exceeding the de Broglie 
wavelength of the excitations. In this case we obtain 

K —f>^ = —wm{(T/n),l} ■ (3-11) 

Thus, we see that the density fluctuations are small at any temperature T <^Td- 

The result of Ea. (|3.11|l looks quite different from what one finds for a uniform ID Bose gas. 
In the latter case, using the Bogoliubov spectrum (|2.33(l and excitations wavefunctions Ij2.32|l . 
for T ^> n we have (8n 2 zz ,)T ~ (T '/ ^/fiT^jn 2 . The difference from Ea. (|3.11|l is related to the fact 
that in the trapped case the density profile shrinks when decreasing temperature well below the 
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global temperature of quantum degeneracy T d = Nftw. Therefore, at T -C Td the local value 
of the degeneracy temperature in the center of the trap is T ~ 2irh 2 n 2 tm /m ~ ^> T d . 

Expressing T d in Ea. (|3.11|) through T we restore qualitatively the uniform-gas result. 

As the low-energy fluctuations of the density are suppressed and high-energy fluctuations 
are also small, for finding the single-particle correlation function we may use the field operator 
in the form (|2.51|l : $(z) = ^Juq{z) exp(i</> p (z)), where the operator <fi z represents the low- 
energy part of the phase fluctuations. Then, similarly to the 2D case, the correlation function 
is expressed through the mean square fluctuations of the phase: 

9l (z, z') = (*V)*(*)) = Vn (*)n (z')exp ; (3 12) 

where 8cf> ZZ ' — 4> P {z) — <fr p (z'). On the basis of Eqs. (|2.49(l and H2.30JI . for the thermal phase 
fluctuations we have 



>)t = 'ff M(j + ^ /2) (Pj(x) - PjWNj . (3.13) 

3 = 1 J 



In contrast to the 2D gas, thermal fluctuations of the phase in ID are mostly provided by the 
contribution of the lowest excitations. A direct calculation of Ea. H3.13|) yields 



4<>T 



4T/i 



In 



(1-ar') (1 + x) 



(1+x') (l-x) 



(3.14) 



For z and z' close to the trap center, the logarithmic term in Ea. H3.14|) is equal to 2\z — 
z'\/Rtf "C 1. Otherwise, it is of order unity, except for the region close to the Thomas-Fermi 
border of the density distribution. 

The temperature at which the quantity {5<p 2 zz ,) is of order unity on a distance scale 
1 2 — z'\ ~ Rtf, is given by 

T^ = T d — . (3.15) 
A* 

Thus, for T <^ both the density and phase fluctuations are suppressed, and there is a true 
condensate. The condition l|3.8|l always provides the ratio T^/fku s» (4iV/a 2 ) 1 / 3 » 1. 

In the temperature range, where T d ^> T ^> T^, the density fluctuations are suppressed, 
but the phase fluctuates on a distance scale 1$ s» Rtf{T<p/T) <C Rtf- Hence, similarly 
to the 2D case, we have a condensate with fluctuating phase (quasicondensate). The phase 
coherence length greatly exceeds the correlation length: 1$ fa l c (T d /T) ^> l c . Therefore, 
the quasicondensate has the same density profile and local correlation properties as the true 
condensate. 

In Fig. EI we present the state diagram for the finite-temperature trapped ID gas at a = 10 
(N* = 100). One clearly sees three quantum degenerate regimes: the BEC regimes of a 
quasicondensate and true condensate, and the regime of a trapped Tonks gas [57] . For N ^> 2V«, 
the decrease of temperature to below T d leads to the appearance of a quasicondensate which at 
T < T$ turns to the true condensate. In the T — N plane the approximate border line between 
the two BEC regimes is determined by the equation (T/hw) = (32N/9N* ) 1 / 3 . For N < 
the system can be regarded as a trapped Tonks gas. 

A cross-over from one regime to another is always smooth. The absence of a sharp transi- 
tion from true to quasiBEC is seen from the behavior of the single-particle correlation function. 
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Fig. 7. Normalized single-particle correlation function gi(z, — z)/no(z) versus \z\/Rtf- 



Using Eqs. (|3.12|) and i|3.14|l and omitting small vacuum fluctuations, for the normalized cor- 
relation function at z' = -z we obtain 

fli (*,-*) = { 1-\z\/Rtf \ 4T/3T * , . 

n (z) \l + \z\/R TF ) ■ V'°> 

This function is displayed in Fig. [7] for various ratios T/T$. In particular, we see that the full 
phase coherence requires temperatures well below T^, 

Present facilities allow one to achieve both the true and quasicondensate regimes for the 
ID trapped Bose gas. For example, in the case of 10 3 rubidium atoms trapped in a cylindrical 
harmonic potential with axial frequency uj ~ 1 Hz, we have the temperature of quantum 
degeneracy Td ~ 100 nK. Then, if the tight radial confinement providing the ID regime is 
made at a frequency ujq ~ 100 Hz, the characteristic temperature is about 10 nK. 

The phase fluctuations lead to a drastic difference of the phase coherence properties of a 
quasicondensate from those of a true condensate. This can be understood from a gedanken 
'juggling' experiment similar to those with 3D condensates at NIST and Munich [58,59]. Small 
clouds of atoms are ejected from the main cloud by stimulated Raman or RF transitions. 
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Observing the interference between two clouds, simultaneously ejected from different parts of 
the sample, allows the reconstruction of the spatial phase coherence properties. Repeatedly 
juggling clouds of a small volume Q from points z and z' of the ID sample, for equal time of 
flight to the detector we have the averaged detection signal / = fl[n (z)+n (z')+2{ , ^'{z)ip{z'))] 
and thus measure directly the single-particle correlation function gi(z, z') = (ip'(z)i/j(z')). 

In the regime of a weakly interacting gas, at T <C T$ the phase fluctuations are small 
and one has a true condensate. In this case, for z' = —z we have (ip^ (z)ip(z')) = rio(z) and 
I = AQhq(z), which shows a pronounced interference effect. The detected signal is twice as 
large as the number of atoms in the ejected clouds. The phase fluctuations grow with T and for 
T > T^, where the true condensate turns into a quasicondensate, the detection signal decreases 
as described by gi(z, —z) from Ea. l|3.I6[l . For T the phase fluctuations completely destroy 

the interference between the two ejected clouds, and / = 2Q,tiq(z). 

The time-dependent single-particle correlation function (w^z, t)^(z', t')) has been calcu- 
lated in [53] and it has been shown how this correlation function can be measured in two-photon 
Rahman outcoupling experiments. 



3.4 Phase coherence in 3D elongated condensates 

The one-dimensional character of thermal phase fluctuations is also present in a very elon- 
gated 3D Bose gas, which leads to the appearance of a quasicondensate in this system. This 
phenomenon, predicted [60] and observed [61-65] in recent studies, is described by the fol- 
lowing physical picture. Excitations of elongated condensates can be divided into two groups: 
"low energy" axial excitations with energies e v < Hljq, and "high energy" excitations with 
e„ > hujQ. The latter have 3D character as their wavelengths are smaller than the radial size 
R±. Therefore, as in ordinary 3D condensates, these excitations can only provide small phase 
fluctuations. The "low-energy" axial excitations have wavelengths larger than R± and exhibit 
a pronounced ID behavior. Just these excitations give the most important contribution to the 
long-wave axial fluctuations of the phase. 

A detailed description of fluctuations in elongated 3D condensates is given in [60], and 
in this Lecture we only briefly outline the results. We consider a cylindrical Thomas- Fermi 
condensate with an axial size Rtf greatly exceeding the radial size R±, and assume that most 
particles are Bose-condensed. The axial thermal fluctuations of the phase, with wavelengths 
larger than R±, are similar to the fluctuations in the ID case. In particular, one finds that the 
mean square fluctuations on a distance scale ~ Rtf are approximately equal to T/T^, where 
the characteristic temperature is given by 

r. = 15 M! W , (3.ir, 

with N being the number of particles, and u> the axial trap frequency. This is qualitatively 
the same as the result of Ea. (|3.15|l where one substitutes Td = Nhus. Accordingly, the phase 
coherence length is again determined by the relation = Rtf^T^/T). It is much larger 
than the healing length: l,f,/l c w (T c /T)(T c /huj ) 2 > I, with T c » (Nufa) 1 / 3 being the BEC 
transition temperature, and ujq the radial confinement frequency (loq 3> ui). 

We thus see that the situation is quite similar to that in ID trapped gases. One has a 
continuous transformation of a quasicondensate into a true BEC when decreasing temperature 
to below T^. 

Most important is the dependence of the phase fluctuations on the aspect ratio of the cloud 
ujq/uj. Fig. IBI shows the ratios T c /T^, /i/T^, and the temperature as functions of loq/us for 
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Fig. 8. The ratios T c /T$ and y,/T$ in (a) and the temperature T$ in (b), versus the aspect ratio uio/oj 
for trapped Rb condensates with N = 10 and too = 500 Hz. 

Thomas-Fermi rubidium condensates at N = 10 5 and ojo = 500 Hz. From these results we see 
that 3D quasicondensates can be obtained in elongated geometries with luq/lu > 50. 

The phase fluctuations are very sensitive to temperature. From Fig. 8 we see that one can 
have T^/Tc < 0.1, and the phase fluctuations are still significant at T < fi, where only a 
tiny indiscernible thermal cloud is present. This suggests a principle for thermometry of 3D 
Bose-condensed gases with indiscernible thermal clouds. If the sample is not an elongated 
quasicondensate by itself, it is first transformed to this state by adiabatically increasing the 
aspect ratio uj p /lo z . This does not change the ratio T/T c as long as the condensate remains 
in the 3D Thomas-Fermi regime. Second, the phase coherence length 1$ or the single-particle 
correlation function are measured. These quantities depend on temperature if the latter is of 
the order of or larger. One thus can measure the ratio T/T c for the initial cloud, which is 
as small as the ratio T^/T c for the elongated cloud. 

Pronounced phase fluctuations have been first observed in Hannover experiments with very 
elongated cylindrical 3D condensates of up to 10 5 rubidium atoms [61,62]. The expanding 
cloud released from the trap was imaged after 25 ms of time of flight and the images showed 
clear modulations of the density (stripes) in the axial direction (see Fig.EJ). The physical reason 
for the appearance of stripes is the following. In a trap the density distribution does not feel 
the presence of the phase fluctuations, since the mean- field interparticle interaction prevents 
the transformation of local velocity fields provided by the phase fluctuations into modulations 
of the density. After switching off the trap, the cloud rapidly expands in the radial direction, 
whereas the axial phase fluctuations remain unaffected. As the mean-field interaction drops to 
almost zero, the axial velocity fields are then converted into the density distribution. 

The mean square modulations of the density in the expanding cloud provide a measure 
of the phase fluctuations in the initial trapped condensate. A direct relation between these 
quantities has been established from analytical and numerical solutions of the Gross-Pitaevskii 
equation for the expanding cloud, with explicitly included initial fluctuations of the phase [61]. 
The obtained phase coherence length was inversely proportional to T, in agreement with theory, 
and for most measurements it was smaller than the axial size of the trapped Thomas-Fermi 
cloud. This implies that the measurements were performed in the regime of quasicondensation. 

The properties of quasicondensates and the phase coherence length were measured directly 
in Bragg spectroscopy experiments with elongated rubidium BECs at Orsay [63] . In this type 
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Fig. 9. Absorption images and corresponding density profiles of BECs after 25 ms time-of-flight in the 
Hannover experiment for aspect ratios [wo/w = 10 (a), 26 (b), 51 (c)]. 



of experiment one measures the momentum distribution of particles in the trapped gas. The 
Orsay studies [63] find a Lorentzian momentum distribution characteristic of quasicondensates 
with axially fluctuating phase [64], whereas a true condensate has a Gaussian distribution. 
The width of the Lorentzian momentum distribution is related to the phase coherence length 
at the trap center. These investigations are described in detail in the lecture of Philippe 
Bouyer. The phase coherence length has also been found in Hannover experiments [65] from 
the measurement of the intensity correlation function of two interfering spatially displaced 
copies of phase-fluctuating condensates. 

It is important to emphasize that the measurement of phase correlations will allow one 
to study the evolution of phase coherence in the course of the formation of a condensate out 
of a non-equilibrium thermal cloud. This problem has a rich physics. For example, recent 
experiments on the formation kinetics of trapped condensates [66] indicate the appearance of 
non-equilibrium quasicondensates slowly evolving towards the equilibrium state. 

Problem: Calculate the momentum distribution for the weakly interacting ID Bose gas in 
a rectangular box of size L at temperatures much smaller than the temperature of quantum 
degeneracy. Describe how the momentum distribution changes when the gas transforms from 
the true to quasicondensate. 
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4 LECTURE 4. CORRELATIONS IN STRONGLY INTERACTING ID BOSE GASES 

In Lectures 2 and 3, we saw that the reduction of spatial dimensionality in trapped Bose 
gases increases the number of possible quantum degenerate states. In 3D we have (true) Bosc- 
Einstcin condensates and they were extensively studied in a large variety of experiments during 
last years. In the 2D case, there are two types of Bose-condensed states: true condensates and 
quasicondensates. In ID, besides these BEC states, one can have a trapped Tonks gas which 
should exhibit some of the fermionic properties. Static and dynamical properties of strongly 
interacting trapped Bose gases are discussed in the lectures of Maxim Olshanii and Sandro 
Stringari. 

In this Lecture we discuss correlation properties of strongly interacting Bose gases, which are 
drastically different from those in the weakly interacting regime. This is especially interesting 
in view of a rapid progress of experimental studies which have already reached an intermediate 
regime between the weakly and strongly interacting ones [67,68]. 

Beyond the weakly interacting regime one can no longer use the mean-field Bogoliubov 
approach. Nevertheless, long-wave properties are generic and can be studied relying on the 
quantum hydrodynamic approach developed by Haldane [48] and widely used for uniform ID 
systems (see [69] for a recent overview). We will employ this approach for studying long- 
distance phase coherence in ID trapped Bose gases [70]. 

The problem of short-range correlations requires more sophisticated approaches. The term 
"short-range" is used in the sense that the distance is of the order of or smaller than the 
characteristic correlation length of the gas, which for repulsive interaction is l c = h/y/mp, with 
\x being the chemical potential. However, the distance is still much larger than the radius 
of interatomic potential, R e . Then, in the purely ID case such short-range correlations may 
be investigated by using the Lieb-Linigcr model which assumes a delta-functional interaction 
between the atoms. 

The Lieb-Liniger model for a uniform system is exactly solvable by using the Bethe Ansatz 
[71]. Thermodynamic functions for this model at zero and finite temperatures have been 
found by Lieb and Liniger [43] and by Yang and Yang [45]. On the other hand, the problem 
of correlation properties is far from being completely resolved, except for some correlation 
functions in the limiting cases of weak and strong interactions. For example, the case of 
infinitely strong interactions is to a certain extent equivalent to that of free fermions and the 
interactions play the role of Pauli principle [42]. In this limit, any correlation function of the 
density is given by the corresponding expression for fermions [72]. The expressions for the 
one-body and two-body correlations for an arbitrary interaction strength were obtained by 
using the Inverse Scattering Method [73] . However, closed analytical results can be found only 
as perturbative expansions in the limiting cases of strong and weak interactions [74-78] . 

We will demonstrate the use of the Lieb-Linigcr model for finding local 2-body and 3-body 
correlations, that is the correlations at distances much smaller than the correlation length 
l c [70,79,80]. These correlations are responsible for the stability of the gas with regard to 
intrinsic inelastic processes, such as 3-body recombination. 

4.1 Lieb-Liniger model for trapped ID gases 

In the case of ID trapped gases, particles undergo zero point oscillations in two (radial) tightly 
confined directions and there is a question to which extent one can use the ID Lieb-Liniger 
model for describing the system. We will consider repulsive interaction between particles. The 
ID regime is realized if the amplitude of radial zero point oscillations Iq = y/ hj mu>o is much 
smaller than the (axial) correlation length l c = fij ^/rnjl, where cjq is the frequency of the 
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radial confinement and fi is the chemical potential of the ID system. This is equivalent to 
the condition /x <C Huiq. One then has a ID system of bosons interacting with each other via 
a short-range potential characterized by an effective coupling constant g > 0. This constant 
is expressed through the 3D scattering length a [56] and for l a ^> a is given by Eq.l j3.lf l. 
Accordingly, the ID interaction (scattering) length is r g ~ l^/a ^> a. In the weakly interacting 
regime, the chemical potential is fi « gn, and the condition l c 3> Iq leads to the inequality 
na < 1. In the Tonks-Girardeau regime the correlation length l c ~ and one should have 
nlo <C 1. We thus see that irrespective of the interaction strength, it is sufficient to satisfy the 
inequalities 

a < l < - . (4.1) 

n 

Then the ID regime is reached and correlation properties of the system can be analyzed on 
the basis of the ID Lieb-Liniger model, which in the absence of an axial trapping potential is 
described by the Hamiltonian 

- fi 2 

ff = 5>2^+s5>(z;-^) , (4-2) 

2=1 i<j 

with N being the number of particles. From the Hamiltonian If4.2f) one easily finds that the 
ratio of the characteristic kinetic energy of particles at the mean interparticle separation to the 
interaction energy per particle, is given by the parameter 7 = mg/h 2 n introduced in Lecture 
3 in Ea. lj3.2fi . In the strongly interacting regime we have 7> 1. 

At zero temperature the energy of the system Eq can be written as 

E h 2 n 2 

/7 = ^r e(7) ' (4 - 3) 

where 7 is given by Ea. (j3.2fi but is not necessarily small, and the function e(7) follows from the 
Bethe Ansatz solution found by Lieb and Liniger [43] . This function is calculated numerically 
for any value of 7 [43,81,82]. In the weakly interacting regime of a quasicondensate (7 <C 1) 
we have 

e(7) = 7 - T^ 3/2 ; 7<1 . ( 4 - 4 ) 

which coincides with the result of the Bogoliubov approach. For the strongly interacting 
Tonks-Girardeau regime (7 ^> 1) the function e("f) is given by 

e W = y(l-^) ; 7»1 ■ (4-5) 
The expression for the chemical potential follows immediately from Eq. lj4.3fl : 



^=w=^i 3e w-^j ■ (4 - 6) 

Accordingly, in the weakly interacting regime (7 <C 1) we have fj, = ng, and in the strongly 
intertacting Tonks-Girardeau regime the chemical potential is 

A 4 = -s — 1 - 5- ; 7 » 1 • (4-7) 



2m \ 37 

In the presence of an (axial) trapping potential Viz) — mui 2 z 2 /2, one should add this 
potential to the Hamiltonian Ij4.2fl . Then, for an arbitrary 7 the model is no longer integrable. 
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However, if the chemical potential of the ID trapped gas satisfies the inequality fio 3> huj, 
one can use the local density approximation [81,83]. Namely, assuming local equilibrium, the 
density distribution is governed by the equation 

fi(n(z)) + V(z) = ti , (4.8) 

where the local value of the chemical potential /j,(n(z)) follows from the solution of the Lieb- 
Linigcr model for a uniform gas with density equal to n(z). In other words, the local chemical 
potential is given by Eg. I|4.6fl in which n = n{z) and 7 = mg/h 2 n(z). 

As jtx(n(z)) should be positive at any local density, Ea. il4.8|) leads to the Thomas-Fermi 
density profile with a (half)size Rtf = (2^0/mw 2 ) 1 / 2 • So, in the interval of distances, —Rtf < 
z < Rtf, the density distribution is governed by Ea. (|4.8|) . and otherwise one has n(z) = 0. 
The density is maximum at the trap center (z = 0) and it smoothly decreases to zero when 
moving to the Thomas-Fermi border of the trapped gas. The normalization condition 




dz n{z) = N 



gives a relation between the chemical potential [io and the number of particles N. 



4.2 Phase coherence at zero temperature 

In Lecture 3 we showed that for a realistic (axial) size of the sample, the ground state of the 
weakly interacting trapped ID Bose gas (7 <C 1 and T = 0) is a true Bose-Einstein condensate. 
The question now is how deeply one should be in the weakly interacting regime that this 
statement is valid. In other words, how small should be the parameter 7 at the trap center for 
having the full phase coherence of the trapped ID gas. 

We will consider a harmonically trapped ID Bose gas in the Thomas-Fermi regime at 
T = and calculate the single-particle correlation function gi(z,z') for distances \z — z'\ 
greatly exceeding the correlation length l c [70]. Assuming a priori small density fluctuations 
at such distances, we will rely on the ID hydrodynamic approach [48] which describes long- 
wave properties of the ID fluid in terms of the density fluctuations 5n and phase 6. Small 
fluctuations of the density at large distances lead to linearized continuity and Euler equations: 

dSn d . . 86 , , . 

at dz az 

^ = -nv N (z)Sn(z) , (4.10) 
where the velocities vj and are given by 

vj = ^ , (4.11) 



N Trh dn 



M • (4 - 12) 

i(z) 



and the ratio K(z) — ^Jv,j(z)/vn(z) is the local value of the Luttinger parameter. Using 
commutation relations l|2.40|l equations (|4.9|l and (|4.10|l can be obtained as equations of motion 
from the quantum hydrodynamic Hamiltonian: 



H = — Id: 

2tt 



(v N (z)(ir5n) 2 + v.j{z){d z & 2 ) ^ , (4.13) 
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where Ej and bj are eigenenergies and annihilation operators of elementary excitations. Note 
that in the uniform case these excitations are phonons and the velocity of sound is c s = ^JvJv~n- 
Similarly, for the trapped gas the local value of the sound velocity is c s — -\/vj(z)vn(z). The 
Hamiltonian (|4.13(l is a generalization of the effective Hamiltonian of Haldane [48] to a non- 
uniform system. 

The solution of Eqs. (|4.9|l and (|4.1U|I is given by the expansion of operators Sh and <fi in 
eigenmodes characterized by an integer quantum number j > 0: 

1 /2 

^ z ,t) = ^-i(^^\ f j (z)b]exp(-is j t/h)+h.c. . (4.15) 
The eigenfunctions are normalized by the condition 



wat (a:) 



where we have introduced a dimensionless coordinate x — z/Rtf- Equations (|4.9|l and (|4.1U|) 
are then reduced to the eigenmode equation: 

^)^-vj(z)^M+ £ l fj (z) = . (4.16) 

Considering fj as a function of the reduced coordinate x — z/Rtf, Ea. (|4.16|l takes the form: 

(1 - a; 2 )// - (2x/p(x))f j + (2e 2 /(3(x)h 2 u; 2 )f 1 = . (4.17) 

The quantity (3{x) — d\n [i/dlnn is determined by the local parameter 7(2) = mg/h 2 n(z). In 
the Tonks regime (7 — * 00) we have (3 — 2, and in the weakly interacting regime (3=1. 

The coordinate dependence of (3 is smooth and we simplify Eq. I|4.17|l by setting (3(z) = 
/3(7o) = A)i where 70 is the value of 7 at maximum density. This simplification has been re- 
cently used [83,84] to study the excitation spectrum of trapped ID Bose gases. Then Eq. (I4.17|) 
yields the spectrum 

e 2 = h 2 Lu 2 (j(3 /2)(j + 2/(3 -l) , 
and the eigenfunctions fj(x) are expressed through Jacobi polynomials: 



f / x _ VU + g + V2)L(, ? + l)L(j + 2oTT) (a , a) 
/JW " 2«r(j+a + l) ■> W ' 

where a = 1/fio — 1. 

As the density fluctuations on a large distance scale are small, for the single-particle corre- 
lation function one can write 

9l (z,z') = = Vn(z)n(z>) exp (~{(4>(z) - ^')) 2 >) . (4.18) 

Using Eq. (|4.15(l the mean square fluctuations of the phase in the exponent of Eq. Ij4.18|) are 
reduced to the sum over j-dependent terms containing eigenfunctions fj and eigenenergies Sj . 
For the vacuum phase fluctuations this sum is logarithmically divergent at large j, which is 
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Fig. 10. Normalized single-particle correlation function gi(z, — z)/n(z) versus z/Rtf, for N = 10 4 
and various values of 70. The solid curves show numerical results, and the dashed curves the results 
of the quasiclassical approach. 



similar to the high-momentum divergence in the uniform case. Accordingly, we introduce a 
cut-off j m ax following from the condition Sj f=a min{/i(z), n(z')} and ensuring a phonon-like 
character of excitations at distances z and z' . This is equivalent to having only the low- 
energy part of the phase fluctuations, ((4> v (z) — 4> p (z')) 2 ) , in the exponent of Ea. (|4.18|l when 
considering the weakly interacting regime. 

The vacuum phase fluctuations have been calculated by using two approaches: numerical 
summation over the eigcnmodcs with exact fj , Sj from the simplified Eq. (|4.17|) , and quasiclas- 
sical approach assuming that the main contribution comes from excitations with j 3> 1. In the 
latter case, for z' — —z we obtain 



which is close to Haldane's result for a uniform system [48] with the Luttinger parameter 
K(z) and correlation length l c (z). In the Tonks regime we have K = 1 and l c = l/n(x) = 
l/n(0)y/l — (x/Rtf) 2 - Remarkably, in this limit the hydrodynamical expression (|4.19() coin- 
cides with the recently obtained exact result [85] for the single-particle correlation function of 
a harmonically trapped Tonks gas. 

The dependence of g\ on the dimensionless coordinate x is governed by two parameters: 70 
and the number of particles N. In Fig.^|we present the quantity gi(z, —z)/n(z) for N = 10 4 
and various values of 70. As expected, the phase coherence is completely lost in the strongly 
interacting regime (70 ^> 1). Moreover, on a distance scale z ~ Rtf the coherence is already 
lost for 70 ~ 1, and the full phase coherence requires 70 well below 0.1. 

4.3 Local correlations at T = 

The strong transverse confinement required for the ID regime can lead to high 3D densities of 
a trapped gas. At a large number of particles the 3D density can exceed 10 15 cm -3 and one 
expects a fast decay due to 3-body recombination. It is then crucial to understand how the 
correlation properties of the gas influence the decay rate. For this purpose, we calculate local 




(4.19) 
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correlations in the ID Bose gas and show that the decay rates are suppressed in the Tonks- 
Girardeau and intermediate regimes, which is promising for achieving these regimes with a 
large number of particles. 

The rate of 3-body recombination is proportional to the local 3-particle correlation func- 
tion g 3 = (^(z)^(z)^(z)^(z)i/(z)^(z)} [86], where all field operators are taken for the 
same time. Similarly, the rates of 2-body inelastic processes involve the correlation function 
52 = (w(z)W(z)^(z)^f(z)). Assuming that local correlation properties are insensitive to the 
geometry of the system we consider a uniform ID gas of N bosons on a ring of circumference 
L. The gas is described by the Lieb-Liniger Hamiltonian 1)4.2(1 which we now rewrite in the 
second quantization: 

h= f dz (^™ + imm) . (4.20) 

J \y2m dz dz 2 J V ' 

For finding g 2 at T = 0, we use the Hellmann-Feynman theorem [88,89]. Namely, one shows 
that the expectation value of the four-operator term in the Hamiltonian (|4.20l) is proportional 
to the derivative of the ground state energy with respect to the coupling constant: 

The first identity follows from the normalization of the ground state wavefunction $o, and 
the second one is obtained straightforwardly from the Hamiltonian ((4.20(1 . The ground state 
energy is given by Ea. ((4.3(1 . and for the 2-particle local correlation, we then obtain 

. 92 (7)=n 2 ^ . (4.21) 

In fact, the original work of Lieb and Linigcr contains a similar calculation of the interaction 
energy. 

The function ^(t)/" 2 calculated by using numerical results for e{^) from [82], is shown in 
Fig- El The quantum Monte Carlo calculations of g 2 [87] arrive at the same results. For small 
and large values of 7, relying on Eqs. 1(4.4(1 and 1(4.5(1 . we obtain: 

^ = i-Vr , 7«i ; (4.22) 

92(l) 4tt 2 

— — = tt~2~ , 7 > 1 . (4.23) 

n z 37^ 

The results in Fig. 1111 and Eq. ((4.23|) clearly show that 2-particle correlations and, hence, 
the rates of pair inelastic processes are suppressed for 7 > 1. This provides a possibility for 
identifying the Tonks- Girardeau and intermediate regimes of a trapped ID Bose gas through 
the measurement of photoassociation in pair interatomic collisions. 

Note that for the weakly interacting regime the result of Ea. ((4.22(l can be obtained directly 
from the Bogoliubov approach. Writing the field operator in the form 1(2.39(1 and using the 
commutation relations ( 12.40(1 . we find 

¥(z)^(z')${z)${z') = h{z)n(z') - hS(z - z') . 

We then represent the operator of the density as n(z) = n + Sn(z) and confine ourselves to the 
second order in small density fluctuations 5h{z) around the mean density n. This gives 

(* t (z)* t (z')*(z)#(z')) =n 2 + {{5h{z)5h{z'))) - n5{z - z') . 
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Fig. 11. Local correlation function g2 versus 7. 



Using the expansion of the density fluctuations in terms of Bogoliubov excitations, given by 
Ea. lj2.46p . we obtain: 

r°° dk 

{&(*)*Hz'M*Mz')) =n 2 + n / — [\f-\ 2 (l + 2N k )-l]exp[ik(z-z')} , 

where E}. = h 2 k 2 /2m, e k = \J E'i + 2fj,E). is the Bogoliubov excitation energy, and the Bogoli- 
ubov function fj~ is determined by Ea. lj2.32p . Then, taking the limit z' — » z we arrive at the 
expression for the local 2-body correlation: 



r°° rlk 

g 2 = (&(z)&(z)*(z)*(z)) = n 2 + nj ^ 



^(l + 2JV fc )-l 



(4.24) 



For T — the occupation numbers for the excitations, Nk — 0, and the integration of Eq. lj4.24p 
immediately gives the result of Ea. (j4.22p . 

The 3-particle local correlation g$ cannot be obtained from the Hellmann-Feynman theorem. 
In the weakly interacting regime (7 <C 1) one can use the Bogoliubov approach, which gives 
Eq. Ij4.24l) with an extra factor 3 in the second term of the rhs. At T = we then have 



53(7) 



(> 



V7 



7 < 1 



(4.25) 



For the Tonks- Girardeau regime (7 3> 1) we demonstrate a method for calculating the 
leading behavior of local correlations. Details are given in [80], and here we present a compact 
derivation of 53 at T = [70]. In the first quantization the expression for this function reads 



53(7) 



Nl 



3\(N- 3)! 



cIza . . .dz 



N 



^ 7) (0,0,0,z 4 , 



zjv, 



(4.26) 



(7) 

where $q is the ground state function given in the domain < Z\ < . . . < z^ < L by the 
Bethc Ansatz solution: 



*o 7) (^i, *2, . ■ ■ , z N ) oc a{P) cxp (i ^ k P 3 



(4.27) 
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where P is a permutation of N numbers, quasimomenta kj are solutions of the Bethe Ansatz 
equations, and 



a(P) = II 

3<l 



i^n — fcp, + k Pl 



For 7 > 1, we extract the leading contribution to $q at three coinciding points by sym- 
metrizing the amplitudes a(P) over the first three elements of the permutation P: 

i^ a (P pi ,P P2 ,P P3 ,P 4 ,...P J v)^-^I]( fc ^-^) > ( 4 - 28 ) 

p K 7 ) }<l 

where j, I = 1, 2, 3. The sign of the permutation P is ep, and p runs over six permutations of 
1, 2, 3. For large 7, the difference of quasi-momenta kj from their values at 7 = 00 is of order 
I/7 and can be neglected. Then, from Eqs. 14.27J1 and (|4.28Jl we conclude that to this level of 
accuracy the ground state wave function at three coinciding points is given by derivatives of 
the wave function of free fermions $q ( z i> z 2, 23, 24, . . .) at z\ = Z2 = 23 = 0: 

$ ( ( 7 ) (0,0,0,z 4 ,...)^-7^3[n(^-^)] <i, OO) ■ (4-29) 

{rin > j<i 

Substituting Eq. H4.29fl into Eq. 1)4. 26[) we express the local correlation g$ through derivatives of 
the 3-body correlation function of free fermions. Using Wick's theorem the latter is given by a 
sum of products of 1-particle fermionic Green functions G(x — y)= J_^ p dk e lk ( x ~ y ) / '2tt , where 
kp = nn is the Fermi wavevector. The calculation from Eq. H4.26fl is then straightforward and 
in the considered limit of 7 > 1 we obtain 



33(7) _ 36 

which yields 



n 3 7677,9 



(G"(0)) 3 -G (4) (0)G"(0)G(0) , (4.30) 
53 16^ 6 



,3 !5 7 6 



(4.31) 



This result for 53 and the result of Ea. H4.23|l for g2, have a transparent physical explanation. 
A characteristic distance related to the interaction between particles is r g — h 2 /mg ~ l/jn, 
and the strong repulsion between particles provides fermionic correlations at interparticle dis- 
tances z > r g . For smaller z the correlation functions practically do not change. Therefore, 
the local correlation g 2 at a finite large 7 is nothing else than the pair correlation function for 
free fermions at a distance r g . The latter is <?2 ~ n 2 (kFr g ) 2 ~ n 2 /^ 2 , which agrees with the 
result of Eq. I|4.23|l . Similarly, is the free-fermion 3-particle correlation function at distances 
~ r g . It is approximately equal to the product of three pair correlation functions, i.e. we have 
33 ~ n e /7 s in qualitative agreement with Ea. H4.23Jl . 

Thus, from Eq. I|4.31|l we conclude that the 3-body decay of ID trapped Bose gases is 
strongly suppressed in the Tonks- Girardeau regime. Moreover, Eq. (|4.25|l shows that even in 
the weakly interacting regime of a quasicondensate, with 7 ps 1CU 2 , one has a 20% reduction 
of the 3-body rate. Thus, one also expects a significant reduction of the 3-body decay in the 
intermediate regime. 

For Iq 3> a, the 3-body recombination process in ID trapped gases occurs at interparticle 
distances much smaller than Iq. Therefore, the equation for the recombination rate is the 
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same as in 3D cylindrical Bose-Einstein condensates with the Gaussian radial density profile 
n 3D = ( n / 7r 'o) ex P ( — /° 2 /'o)i where p is the radial coordinate. There is only an extra reduction 
by a factor of g 3 /n 3 . A characteristic decay time r is then given by the relation 



n fo j f n 3D\ 3 a 3D g 3 

- = 93 a 3D J 2npdp (—) ^^-^ 



(4.32) 



where ol 3 d is the recombination rate constant for a 3D condensate. Even for 3D densities 
u/ttIq ~ 10 15 cm -3 the life-time r can greatly exceed seconds when approaching the Tonks- 
Girardeau regime. 

In the recent rubidium experiment at NIST [68], an array of ID tubes of bosons was created 
by optically confining the atoms with a radial frequency ujq ~ 30 kHz (Iq w 350 A). Then, for 
the number of atoms N s» 200 in each ID tube, the intermediate regime with 7 ~ 1 has 
been achieved. The measurement of the 3-body decay showed a reduction of the rate by 
approximately a factor of seven. 



4.4 Finite-temperature local correlations 

We now discuss local correlations in ID Bose gases at finite temperatures and confine our- 
selves to two-body correlations in the uniform case [79]. The calculation of g% will allow us to 
identify the various regimes of quantum degeneracy, which is important for the development of 
possible atom lasers in ID waveguides. At zero temperature there are two physically distinct 
regimes of quantum degeneracy: the weakly interacting regime of a quasicondensate for 7 <C 1, 
and the strongly interacting Tonks-Girardeau regime where the ID gas undergoes "fermion- 
ization" . These regimes are also present at finite temperatures. However, we will see that 
for a very small interaction strength one can have a decoherent quantum regime, where the 
fluctuations are enhanced and reach the non-interacting Bose gas level with g 2 — ► 2n 2 (rather 
than 52 — ► n 2 characteristic for the quasicondensate regime). We will also see that for a large 
interaction strength the reduction of local correlations is present at temperatures exceeding the 
temperature of quantum degeneracy T4 — 2n h 2 n 2 /m and, in this respect, one has the regime 
of " high-temperature fermionization" . 

So, we again consider a ID uniform gas of N bosons on a ring of circumference L, described 
by the Hamiltonian 1)4. 2|) or, equivalently, by the Hamiltonian 1|4.20[) . The exact solution of this 
problem at finite T has been found by Yang and Yang [45] using the Bcthe Ansatz. They derived 
exact integral equations for thermodynamic functions and proved their analyticity. Note that in 
trapped gases at finite temperatures, the ID regime requires the thermal de Broglie wavelength 
of particles Ay to be much smaller than the amplitude of zero point oscillations in the tightly 
confined directions, Iq. This requirement should be added to Ea. (|4.1|) . and we thus obtain 
the inequalities a <C Iq <C Ax}. They allow us to analyze finite-temperature correlation 

properties of the ID trapped gas at any 7 by using the Yang- Yang results. 

The behavior of the 2-body local correlation at finite temperatures is governed by two 
parameters: 7 and r = AirT/Td. We first show how gi has been calculated by using the 
Hellmann-Feynman theorem [88,89]. Consider the partition function Z = exp(— F/T) = 
Trexp(— H/T) which determines the free energy F. Here the trace is taken over the states 
of the system with a fixed number of particles in the canonical formalism or, for the grand 
canonical description, one has to replace the condition of a constant particle number by the 
condition of a constant chemical potential p and add the term —pN to the Hamiltonian. For 
the derivative of the free energy with respect to the coupling constant one has 
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Fig. 12. Correlation function g 



(2) _ 



52 /n 2 versus 7 at different r = 4irT/Td. 



Introducing the free energy per particle / (7, r) = F/N , the 2-particle correlation function is: 



In Fig. 12 we present the results for gi, found by numerically calculating /(7,r) from the 
Yang- Yang exact integral equations and then using Ea. (|4.34|l [79]. 

We now give a physical description of different regimes determined by the values of the 
coupling constant 7 and the reduced temperature r [79]. As well as at T = 0, in the strongly 
interacting regime the correlation function g% reduces dramatically due to the strong repulsion 
between particles, and the physics resembles that of free fcrmions. Interestingly, this regime 
is realized both below and above the quantum degeneracy temperature. Similarly to (73 and 
gi at zero temperature, the finite-temperature 172 can be expressed through derivatives of the 
Green function of free fermions G(x) = J dknp(k) exp (ikx)/(2n), where np(k) are occupation 
numbers for free fermions. For the correlation function we obtain 



In the regime of quantum degeneracy, t -C 1, the local correlation is dominated by the 
ground state distribution np(k) = 6(kp — fc 2 ), where kp — irn is the Fermi momentum, and 
one only obtains a small finite-temperature correction to the zero-temperature result l|4.23f) . 

In the temperature interval 1 <C t -C 7 2 the gas is non-degenerate, but the interaction 
length r g is still much smaller than the thermal de Broglie wavelength At- As well as at 
zero temperature, 32 can be viewed as the pair correlation function for free fermions at a 
distance ~ r g . Taking into account that the characteristic momentum of particles is now the 
thermal momentum kp — y2mT/h, one finds ~ n 2 (fcrr g ) 2 ~ n 2 T/j 2 . Calculating the 
Green function G(x) for the classical distribution np(k), we obtain: 




(4.34) 



4 r 



(G'(0)) 2 -G"(0)G(0) 




K r < 7 2 



(4.35) 
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which agrees with the given qualitative estimate. The correlation function g2/n 2 is still much 
smaller than unity and we thus have a regime of "high-temperature fermionization" . The result 
of Eq. i|4.35|) agrees with the outcome of numerical calculations. 

In the weakly interacting regime of a quasicondensate, one may use Ea. (|4.24|l . For tem- 
peratures t>7, thermal fluctuations are more important than vacuum fluctuations, and we 
obtain: 

^§ = 1 + ^ , 7 «T«V7«1 • (4-36) 

The phase coherence length is determined by thermal long- wavelength fluctuations of the phase. 
The calculation, similar to that for a trapped gas in Lecture 3, gives « H 2 n/mT. The condi- 
tion l(j,^> l c , which is necessary for the existence of a quasicondensate and for the applicability 
of the Bogoliubov approach, immediately yields the inequality r <C ^7. Thus, Eq. 14.36JI is 
valid under the condition 7<r< ypy, and the second term on the rhs of this equation is a 
small correction. In the region of its validity, the result of Eq. 1)4.36(1 agrees well with numerical 
data. 

At a very weak coupling strength given by 7 < r 2 < 1, the gas is in a decoherent quantum 
regime [51,79], where both phase and density fluctuations are large and the local correlation is 
always close to the result for free bosons, 32 = 2n 2 . The only consequence of quantum degener- 
acy is the quantum Bose distribution for occupation numbers of particles, and the decoherent 
quantum regime can be treated asymptotically by employing a standard perturbation theory 
with regard to the coupling constant g. 

The presence of the quantum decoherent regime helps to explain the apparent contradiction 
that in the thermodynamic limit a free Bose gas at any finite temperature is known to display 
large thermal (Gaussian) density fluctuations with 172 = 2rt 2 . For the 3D gas this result 
requires the grand canonical description [90], whereas in ID and 2D it is valid for any choice 
of the ensemble. On the other hand, the widely used mean-field Bogoliubov approach for 
an interacting Bose gas leads to 32 ~ n 2 . The above discussed results explain this fact in 
the exactly solvable ID case: there is a continuous transition from the quasicondensate to 
decoherent regime, depending on the density and temperature. As 7 is decreased towards a 
free gas, the quasicondensate result of <?2 — n 2 only holds above a certain interaction strength. 
Below this, there is a dramatic increase in fluctuations, with 172 —> 2n 2 in the free gas limit. 

Problem: Prove that the density fluctuations for the uniform Tonks-Girardeau gas are small 
in the limit of large distances. Find the distance scale on which they become large. 
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